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Abstract 

The history of the discovery of bosonic string theory is well documented. 
This theory evolved as an attempt to find a multidimensional analogue of Eu- 
ler's beta function to describe the multiparticle Veneziano amplitudes. Such an 
analogue had in fact been known in mathematics at least in 1922. Its mathemat- 
ical meaning was studied subsequently from different angles by mathematicians 
such as Selberg, Weil and Deligne among others. The mathematical interpre- 
tation of this multidimensional beta function that was developed subsequently 
is markedly different from that described in physics literature. This work aims 
to bridge the gap between the mathematical and physical treatments. Using 
some results of recent publications (e.g. J.Geom.Phys.38 (2001) 81 ; ibid 43 
(2002) 45) new topological, algebro-geometric, number-theoretic and combina- 
torial treatment of the multiparticle Veneziano amplitudes is developed. As 
a result, an entirely new physical meaning of these amplitudes is emerging: 
they are periods of differential forms associated with homology cycles on Fer- 
mat (hyper)surfaces. Such (hyper) surfaces are considered as complex projective 
varieties of Hodge type. Although the computational formalism developed in 
this work resembles that used in mirror symmetry calculations, many additional 
results from mathematics are used along with their suitable physical interpre- 
tation. For instance, the Hodge spectrum of the Fermat (hyper) surfaces is in 
one-to-one correspondence with the possible spectrum of particle masses. The 
formalism also allows us to obtain correlation functions of both conformal field 
theory and particle physics using the same type of the Picard-Fuchs equations 
whose solutions are being interpreted in terms of periods. 
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1 Introduction 



1.1 Some properties of Veneziano and Virasoro ampli- 
tudes 

In 1968 Veneziano [1] had postulated 4-particle scattering amplitude A(s,t,u) 
given (up to a common constant factor) by 

A{s, t, u) = V{s, t) + V(s, u) + V(t, u) (1.1) 

where 

l 

V(s,t) = J x-^-^l - x^^dx = B(-a{s),-a{t)) (1.2) 
o 

is Euler's beta function and a(x) is the Regge trajectory usually written as 
a(x) = a(0) + a'x with a(0) and a' being the Regge slope and the intercept, 
respectively. In case of Lorentzian metric with signature {+, — , — , — } the Man- 
dclstam variables s, t and u entering the Regge trajectory are given by [2] 

8 = ( P l+qi)\ (1.3) 
t = {Pl-V2) 2 , 
u = {pi - q 2 ) 2 . 

The 4-momenta pi and qi are constrained by the energy-momentum conserva- 
tion pi + qi = P2 + 92 leading to relation between the Mandelstam variables: 

4 

s + t + u = 'Y^m 2 . (1.4) 
i=i 

Already Veneziano [1] had noticed 2 that to fit the experimental data the Regge 
trajectories should obey the constraint 

a(s) +a(t) +a(u) = -1 (1.5) 

consistent with Eq.(1.4). He also noticed that the amplitude A(s,t,u) can be 
cquivalcntly rewritten with help of this constraint as follows 

A(s,t,u) = r(-a(s))r(-a(t))r(-a(u))[sin7r(-Q!(s))+sin7r(-a(i))+sin7r(-a(u))]. 

(1.6) 

This result looks strikingly similar to that suggested a bit later by Virasoro [3] . 
Up to a constant it is given by 

_ r(a)r(6)r(c) 

A(s '^ u >~ T (a + b)T(b + c)T(c + a) {1J > 



2 To get our Eq.(1.5) from Eq.7 of Veneziano paper, it is sufficient to notice that his 1 — a(s) 
corresponds to ours -a(s). 
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with a = — |a(s),etc. also subjected to the constraint: 



Use of the formulas 



and 



^(a(s)+a{t)+a(u)) = -l. (1.8) 



r(x)r(i - x) = (i.9) 

sin 7tx 



4 sin x sin y sin z = sin(x + y — z) + sin(y + z — x) + sin(z + x — y) — sin(a; + y + z) 

(1.10) 

permits us to rewrite Eq.(1.7) in the alternative form (up to unimportant con- 
stant): 

A(s,t,u) = [T(-^a(s))T(-^a(t))T(-^a(u))} 2 x (1.11) 

[sin7r(-^a(s)) + sin7r(-^a(t)) + sin7r(-^a(u))]. 

Although these two amplitudes look deceptively similar, mathematically, they 
are markedly different. Indeed, by using Eq.(1.6) conveniently rewritten as 

A{a,b,c) =r(a)r(fr)r(c)[sin7ra + sin7rfr + sin7rc] (1.12) 

and exploiting the identity 

wz n z 1 C(l — z ) / , „n 

cos T^rlV (L13) 

after some trigonometric calculations the following result is obtained: 

A(a,b,c) = — m (1.14) 

provided that 

a + b + c=l. (1.15) 

For the Virasoro amplitude, apparently, no result like Eq.(1.14) can be obtained. 
The differences between the Veneziano and the Virasoro amplitudes are much 
more profound as the rest of this paper demonstrates. 

The result, Eq.(1.14), is remarkable in the sense that it allows us to interpret 
the Veneziano amplitude from the point of view of number theory, the theory 
of dynamical systems, etc. Such interpretation is presented in some detail in 
Section 3.1. while Section 2 provides the necessary mathematical background. 

Moreover, in our previous work, Rcf. [4], the following partition function 
describing (train track) dynamics of 2+1 gravity was obtained 

m = ^WT' (L16a) 
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This function is meant to describe a dynamical transition from the pseudo 
Anosov (Appendix B) to the Seifcrt fibered dynamical regime controlled by 
the temperature-like parameter (3. The exact physical nature of the parameter 
(3 was left unexplained. The train tracks " live" on the surface of the punctured 
torus-the simplest surface of negative Euler characteristic on which they can 
"live" [5]. Normally, torus without puncture (s) can be foliated by vector fields 
without singularities. The punctured torus (via Schottky double construction) 
can be converted into a double torus, that is into the Riemann surface of genus 
2. If one tries to foliate such surface with some line field (not to be confused 
with the vector field [6]), one can easily discover that singularities of such field 
will occur inevitably. The number of singularities is controlled by the Poincare- 
Hopf index theorem 3 . Foliations with singularities are known as pseudo Anosov 
(Appendix B) type. For the punctured torus only 2 singularities are allowed and 
their scattering on each other is very much like that which the Veneziano ampli- 
tude is describing. The Mandelstam variables s,t,u can be used in the present 
case as well thus replacing our dynamical zeta function, Eq.(1.16a), with the 
product of zeta functions, i.e. with the Veneziano amplitude, Eq.(1.14). Such 
replacement is not complete however since our dynamical zeta function has 
"wrong" argument: (3 — 1, in the numerator of Eq. (1.16a), instead of 1 — (3. 
Fortunately, this deficiency is easy to correct if we recall the presentation of 
the Riemann zeta function as infinite product over primes. Then, instead of 
Eq.(1.16a), we obtain 

where the product is taken over all prime numbers. Such presentation invokes 
immediately connections with the p-adic string theory to be described below. 
For now, however, we still need to make several observations. To begin, let us 
introduce the p-adic analogue of zeta function, in particular, 

to be reobtained in Section 2, Eq.(2.15). We can equivalently present it as 

W !) = ^T-^T = (-V-Xa - fl- (1-17) 

It should be clear from Eq.(1.17) that Q p {(3 — 1) has the same pole as Q p {\ — (3) 
and, hence, the Veneziano amplitude, Eq.(1.14)> * s obtainable from dynamics of 
2+1 gravity and, therefore, from the full Einsteinian gravity. Such arguments, 
although plausible, not quite rigorous yet. Indeed, in the case of the open string 
the world sheet is the Poincare upper half plane (or, equivalently, the open 
disc). The incoming/outgoing particles are represented as "insertions on the 

3 It is possible to have different sets of singularities as long as their total index remains the 
same [5,6] 
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boundary" in terminology of Ref.[7], page 46, or as cusps in mathematical ter- 
minology. In the case of 4-particle Veneziano amplitude there are 4 "insertions 
on the boundary" . But, as Fig. 3 of our work on dynamics of train tracks, Ref. [4] , 
indicates, in the case of the punctured torus there arc in the simplest case also 4 
insertions at the boundary. Moreover, the open disc Poincare model is not the 
model of the world sheet in the case of 2+1 gravity. It is used for the description 
of the dynamics of simplest train tracks. This dynamics is taking place in the 
Teichmiiller space of the punctured torus which is just the open Poincare disc 
[8] . It can be shown [9] , that the Teichmiiller space of the punctured torus is the 
same as that of the 4 times punctured sphere (considered to be the world sheet 
for the closed string theory, Ref. [7], page 35). Hence, already at the tree level in 
string theory the description of physical processes actually takes place not in real 
but in the Teichmiiller space. And, of course, mathematically rigorous study 
of the dynamics of train tracks indicates that, indeed, such dynamics is taking 
place in the Teichmiiller space. The above picture, perhaps, can be generalized 
to surfaces of genus higher than one 4 . In this paper, however, no attempts are 
made to do this in view of the fact that there is much more mathematically 
elegant and efficient way to accomplish the task. Before discussing details of 
this other way, we would like to make several comments regarding connections 
of the result, Eq.(1.14), with the p-adic string theory. 

1.2 Connections with the p-adic string theory 

Eq.(1.14) was used implicitly (but essentially) in the p-adic string theory as can 
be seen from the review by Bekke and Freund, Ref. [11], and is discussed below. 
Since Eq.(1.14) was also obtained in our earlier work, Ref. [4], in connection 
with dynamics of train tracks representing dynamics of 2+1 gravity, from the 
standpoint of an external observer, the evolution of the train track patterns 
takes place in normal 2+1 or even 3+1 physical space-time 5 . Eq.(1.14), when 
it is obtained with help of the p-adic bosonic string formalism, " lives" instead in 
26 dimensional space- time ( page 26 of Ref. [11]). Hence, it is worth discussing 
and comparing in some detail both cases now. 

Let us begin with the case of the p-adic string theory. It is commonly believed 
that the only rationale of such theory lies in regularizing of the string world 
sheet by placing such sheet on the Bruhat-Tits (known in physics literature as 
Bethc lattice) tree (e.g see Appendix C) "which can be embedded in ordinary 
real space-time" (page 29 of Ref. [11]). Consistency with known continuous 
results in bosonic string theory require, however, that such real space-time is 26 

And, indeed, in a series of papers by Takhtajan et al (e.g see [10] for the latest paper 
and references therein) a serious attempt has been made to accomplish this task. To our 
knowledge, the multiparticle Veneziano amplitudes have not been reobtained thus far in these 
papers. 

5 Moreover, such picture is not just imaginary. It can be made entirely real if one looks at 
dynamics of disclinations in 2 dimensional liquid crystals. Such liquid crystal dynamics is in 
one to one correspondence with the dynamics of train tracks [5,6]. 
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dimensional. Such conclusion is reached roughly speaking based on the following 
arguments. Consider the identity 

r( S -) C(s) = n-^T(^) C(l - a). (1.18) 

Following Rcf.[ll] we define the adelic zeta function Ca( s ) by 

( A (s)=n-ir( S -)((s) (1.19) 

then, Eq.(1.18) can be rewritten as follows 

Ca(«) = Ca(1 - «)■ (1-20) 
Using the same reference, an adelic gamma function Ta(s) is then defined by 

r A (s)= c f {s) (i.2i) 

Ca(! - s ) 

In view of Eq.(1.20), clearly, Ta(s) = 1. This result is a particular example of 
the so called "product formula" well known in the algebraic number theory [12]. 
Following general rules of p-adic analysis beautifully explained in the classical 
book by Artin [12], one can rewrite the statement Ta(s) = 1 in a form of the 
product over both finite and infinite places, i.e. 

r A («) = n r p(«) = i> ( L 22) 

p 

where p runs over all primes (just like in the case of Riemann zeta function) 
and, in addition, over the "prime at infinity". In our case 



1-P- S ~ Cp(l-s) 



but 

Toois) = /~ (S) , = n- s+ ir( S -)/T(^), (1.24) 
Coo(!- s ) 2 2 

where Coo( s ) = ^ ^T(|) in view of Eq. (1.18). Clearly, we get then 

rfr 1 ^ n r P w = L ( L25 ) 



This leads to 



,^c^ = Ii r ' (s)= ^>r = «^> (L26) 
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in view of Eq.(1.18). But the l.h.s. of Eq.(1.26) is obtainable from the p-adic 
string path integral [11, 13] 6 . Hence, the p-adic open string scattering amplitude 
is just an inverse of the standard Veneziano amplitude in view of Eq.(1.14). 
Since the standard amplitude is obtainable via Polyakov path integral in 26 
dimensions, the p-adic string should also live in 26 dimensions in view of the 
product formula, Eq.(1.22). This is the conclusion reached in Ref.[ll]. For- 
tunately (or unfortunately!), this conclusion is incorrect. Indeed, let us take 
another look at Eq-s (1.20)-(1.22). We have 

ncp( s )=n^( i - s ) ( l2? ) 

p p 

or, equivalently, 

CooM n w = Coc(i - s) n c P (i - (i.28) 

p^oo p^oo 

Evidently, it is only natural to expect that "locally" we should also have 

C P («) = C P (1 - s) (1.29) 

for all pis, including p=oo. This is not the case, however, in view of definitions 
of local zeta functions, e.g. see Eq.(1.23). This negative observation cannot be 
simply repaired. Indeed, from p-adic analysis [14] and, even more advanced 
theory (which includes p-adic analysis as its part) of Drinfeld modules [15], it 
can be found that the existing in mathematics (e.g. see Eq.(3.45) below) expres- 
sions for the p-adic gamma functions (actually, there are several expressions 
[15]) differ markedly from that given in Eq.(1.23). Moreover, the Bruhat-Tits 
tree used for calculations of Veneziano and Virasoro amplitudes in Ref-s[ll,13] 
is the discrete analogue of the upper Poincare' half plane [ 16, 17] which is the 
universal covering space for any hyperbolic Riemann surface while the con- 
ventional string theory formulations producing at the tree level Veneziano and 
Virasoro-Shapiro-likc amplitudes involve fiat metrics which is obviously nonhy- 
perbolic. Hence, the whole chain of arguments leading to the conclusion that 
p-adic strings are "living" in 26 dimensions is apparently incorrect. Moreover, 
we argue in this work that use of the Bruhat-Tits tree in the p-adic string the- 
ory is motivated not by the necessity of regularization of the world sheet. It 
occurs rather naturally and is very closely associated with complex multipli- 
cation theory developed by Shimura and Taniyama [18]. It can be thought as 
some analogue of the Fourier component in the Fourier expansion 7 of some con- 
tinuous function. Just like in the linear Fourier analysis our knowledge of one 
Fourier component is practically sufficient for restoration of the whole function, 
the knowledge of one p-adic component is sufficient for restoration of the whole 

6 Actually, the product of 3 factors corresponding respectively to s, t and u channels is 
obtained. 

7 More exactly, it occurs as an analogue of the individual term in Euler's prime number 
product decomposition of zeta function. 
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function as well. This statement is illustrated below on examples of calculation 
of the p-adic analogue of the Veneziano amplitude discussed in Sections 3.3. ,4.3. 
and 5.3.1. 

1.3 Organization of the rest of this paper 

The history of development of the dual resonance models leading ultimately 
to models of relativistic bosonic and fermionic strings is well documented and 
can be found, for example, in the excellent collection of review articles edited 
by Jacob, Ref.[19]. This development one way or another stems from earlier 
efforts to develop the axiomatic S-matrix theory-the idea which can be traced 
back to two papers by Heisenberg [20] of 1943. The mathematical foundations of 
the S-matrix theory can be traced back to even much earlier works of Kramers 
and Kronig dating back to years 1926 and 1927. The usefulness of their results 
(related originally to the light scattering) to particle physics scattering is well 
documented in the classical monograph by Bjorken and Drell, Ref.[21]. The 
contribution of Kramers and Kronig to scattering problems lies in their obser- 
vation of the usefulness of the Cauchy integral formula to such problems. Surely, 
from the modern perspective, one can talk cquivalcntly about usefulness of the 
Ricmann-Hilbcrt problem [22] to scattering problems. And this observation 
leads us ultimately to the Picard-Fuchs type equations considered in Section 
5 . If one reformulates the Kramers-Kronig results into the language of particle 
physics as it is done by Bjorken and Drell [21], then in few words the essence of 
the scattering problem can be formulated as the statement about the scattering 
amplitude 

Im/(w) = ^-<7(w) (1.30) 

connecting the imaginary part of the forward scattering amplitude Im/(o>) with 
the total crossection a{ui) for the "frequency" u> > 0. Suppose, we know the 
amplitude f(to), then, the familiar trick (which is just a corollary of the Cauchy 
theorem) tells us that 

f(u + is) - f(uj - is) = 27nlm/(w) , e -> 0+ (1.31) 

thus allowing us to determine the total crossection using Eq.(1.30). Alterna- 
tively, if such a crossection can be obtained experimentally, one can restore 
the real part of the amplitude through Kramers-Kronig relations thus recon- 
structing the whole amplitude from the experimental data. Such simple minded 
picture becomes very complex when, for example, one is trying to analyze the 
analytic properties of the multiparticle vertex parts (contributing to scattering 
amplitude) based on analysis of the appropriate Feynman diagrams of quan- 
tum field theory. At the physical level of rigor this problem was studied by 
Landau [24]. Subsequently, much more sophisticated cohomological analysis of 

8 The connection between the Picard-Fuchs equations and the Riemann-Hilbcrt problem is 
discussed, for example, in Rcf. [23] 
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the whole problem was developed by Pham [25] whose results were brought to 
perfection by Milnor[26] and Brieskorn [27]. In physics literature the results of 
Pham were carefully analyzed in two fundamental papers by Ponzano, Regge, 
Speer and Westwater [28]. Unfortunately, subsequent development of particle 
physics went into different direction(s). Nevertheless, the major problem of re- 
construction of the underlying particle physics model from scattering amplitudes 
with known analytic properties had ultimately created what is known today as 
string/membrane theory. As it was mentioned in Section 1.1., Veneziano had 
not derived his 4-particle amplitude theoretically. This amplitude was postu- 
lated and subsequently checked experimentally [7,19]. The attempts to extend 
this amplitude to multiparticle scattering had lead ultimately to formulation of 
string theory. Since mathematically, the simplest Veneziano amplitude is just 
Euler's beta function, the extension of such amplitude to the multiparticle case 
mathematically is reduced to the correct multidimensional extension of Euler's 
beta function integral. Such an extension can be found in the book by Edwards, 
e.g. see Ref [29], page 167, published inl922.His result was rediscovered by many 
mathematicians. The latest in this list of people, perhaps, Varchenko whose two 
papers [30] of 1989-90 contain the result of Edwards and, surely, much more. 
Ironically, in 1967-a year earlier than Veneziano paper was published, the paper 
by Chowla and Selberg had appeared [31] which related Euler's beta function to 
the periods of elliptic integrals. The usefulness of results of Chowla and Selberg 
is discussed in Section 4 within the context of conformal and particle physics the- 
ory. The result of Cowla and Selberg was generalized by Andre Weil whose two 
influential papers on the same subject [32, 33] had brought into picture periods 
of the Abelian varieties, Hodge rings, etc. thus inspiring Benedict Gross to red- 
erive Edwards result in 1978 [34] without actually being aware of its existence. 
In the paper by Gross the beta function appears as period of the differential 
form on the Jacobian of the Fermat curve. His results as well as those of Rohlich 
(placed in the appendix to Gross paper) have been subsequently documented 
in the book by Lang [35] . These results are explained from physical standpoint 
in Section 3.2. where they are reinterpreted with help of Milnor's fundamental 
results on singular points of complex hypersurfaces [26]. Milnor's work to a 
large extent had shaped up the theory of singularities and, not surprisingly, all 
the results of this paper could be considered as some practical applications of 
the singularity theory as described, for example, in the monograph by Arnol'd, 
Varchenko and Gussein-Zade [36]. 

Although in the paper by Gross the multidimensional extension of beta func- 
tion is considered, e.g. see page 207 of Gross paper [34], the details were not 
provided however. The details are provided (to our knowledge) only in this 
paper. They are based on ideas developed in lectures by Deligne [37] delivered 
in 1978-1979. Deligne was seemingly unaware of both results of Edwards and 
Gross. In developing of Hodge theory Deligne noticed that Hodge theory re- 
quires some essential changes(e.g. introduction of mixed Hodge structures, etc.) 
in the case if the Hodge-Kahler manifolds possess singularities. We discuss 
these modifications briefly in Section 5 and in Appendix D within the context 
of standard singularity theory. Although this paper is an attempt to present 
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a balanced treatment of both the number-theoretic (including the p-adic ) and 
standard analytic aspects of the problem of periods associated with Veneziano 
amplitudes, and in spite of the fact that the paper came out rather long, in 
reality, it presents just a gentle scratch on the surface of enormous amount of 
information. For instance, the number-theoretic aspects of the mixed Hodge 
structures involving theories of motives [37,38], of crystalline cohomology [37- 
39], of Tannakian categories [37-38], etc. are left completely outside the scope of 
this paper mainly because their physical meaning still remains to be uncovered. 
Also, the symplectic aspects of our scattering problem [40] needed for develop- 
ment of correct "string" -theoretic model associated with Veneziano amplitudes 
are left outside the scope of this paper as well. The emerging opportunity to " fill 
in the gaps" and to demonstrate that the correct formulation of "string" theory 
leads indeed to " the theory of everything" , hopefully, should provide an inspi- 
ration for our readers. In Section 6 we provide some outline of what lies ahead. 
Very exciting connections between many branches of so far " pure" mathematics 
and "down to earth" and, hence," dirty" physics are yet to be unveiled. 

2 Number fields versus function fields 
2.1 Valuations 

The existing physical literature uses real, complex and p-adic numbers. Use of p- 
adic numbers is still rather exotic [11,13, 41] and to a large extent artificial and, 
hence, unpopular. It is strange that function fields have not been used in physics 
so far. In algebra [12,42] and number theory [12,43] it was recognized long ago 
that function fields have practically the same properties as number fields but in 
many cases are easier to handle. We would like to argue that the function fields 
are not only easier to handle than the number fields but also that these fields are 
physically more relevant. Not surprisingly, these fields found their widespread 
use in discrete mathematics, computer science and coding theory [44.45]. 

Elementary number theory begins with the ring of integers Z . This ring 
can be made into the field of fractions Q if the division operation is included 
(provided that there is no divisions by zero) . The field of rational numbers R 
is obtained by the completion operation involving the Cauchy sequences and 
the usual absolute values [46]. The field of complex numbers C is obtained by 
algebraic extension of Z, Q or R. E.g. solution of equation x 2 + 1 = is the 
simplest basic example leading to complex numbers. Evidently every number n 
€ Z can be decomposed into primes by the property of unique factorization, e.g. 

n = ±pr--- P T r , (2.1) 

etc. At the same time, every polynomial f(T) can also be uniquely decomposed 
over the algebraically closed field C as follows: 

f(T) = c(T - ai) mi • • • (T — a r ) m r . (2.2) 
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Comparison between Eq.s (2.1) and (2.2) already suggests a close analogy. It 
is not our purpose here to reproduce long excerpts from algebra related to 
polynomial rings F[T] over the field F where F can be either one of the fields 
above or finite field of q elements. In the last case one usually writes F g [T]. 
Evidently, the function ring F[T] is analogous to the number ring Z while the 
quotients of polynomials taken from the ring F[T] form the field F(T) analogous 
to Q,etc. Computationally the function field has many advantages over the 
number field. We shall discuss only those which are related to the goals we have 
in mind. To this purpose, following Artin [12], we introduce the following 

Definition 2.1. A valuation of a field F is a real-valued function \x\, 
defined for all x E F satisfying the following requirements: 

a) \x\ > 0; \x\ = if and only if x = 0; 

b) \xy\=\x\ \y\ ; 

c) \x + y\ < \x\ + \y\ . 

There are many ways to introduce valuations explicitly as explained in 
Rcf.[47] and we shall take advantage of this fact. Nevertheless, every valuation 
is equivalent to a valuation for which the triangular inequality holds. The trian- 
gular inequality can be sharpened for the non-archimedean valuations. In this 
case instead of c) one has the following inequality 

\x + y\ < max(|a;| ,|y|). (2.3) 

The valuation is called trivial if \x\ = 1 \/x E F* (F* = F \ 0).ln the case of field 
Q it is convenient to choose some prime number p so that any number a E Q, 
a ^ 0,can be presented in the form a = p v b with numerator and denominator 
of b being prime to p. 

Definition 2.2. The standard valuation at p is defined according to the 
rule: 

\a\ p = p~'\ where v — ord p (a) = order of a at p. (2-4) 

Clearly, with such definition the whole field Q is divided into equivalence classes. 
Within such scheme the usual absolute value |o| = H Poo - This rule serves as 
definition of a fictitious " prime at infinity" . The difference between the prime 
at infinity and the rest of primes is profound and represents the major difficulty 
in number theory as explained by Mazur, Ref. [48], page 15. 

Fortunately, for the case of function fields this problem does not exist alto- 
gether since the prime at infinity is treated as any other prime. Let F[T] be a 
polynomial ring and F(T) its quotient field of rational functions. Let p = p(T) 
be an irreducible polynomial (of degree > 1) with leading coefficient 1 in F[T]. 
For each rational function ip(t) E F(T) one can write 

v (T) = {p{T)y f -^± (2.5) 

where f(T) and g{T) E F[T] are polynomials and p \ fg. By analogy with 
Eq.(2.4) the standard valuation can be defined now as 

I^L = cxp{-^dcgp}. (2.6) 
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In this case the analogue of the absolute value (prime at infinity) is just 

Moo =exp{deg(^}. (2.7) 

These results can be made physically relevant if one considers the field of mero- 
morphic functions on the Riemann sphere, i.e. on C+oo. Let f(x) be such 
function defined on C. Then, near some point xq it can be presented as 

f(x) = (x-x ) ord *° f{x) g(x), (2.8) 

provided that g(xo) ^ or oo. Using these results one can write 

\f( x )\ X0 = exp{-z/} where v = ord Xo f(x). (2.9) 

If x =oo, one writes f(x) — (I) ordoo -f( x ) g(I) jW here g(0) ^ or oo, and using 
Eq.(2.7) one defines valuation at infinity as 

|/( a; )| 00 =ex P {deg/( a; )} (2.10) 

since in this case ord 00 f(x) = deg/(x). For the meromorphic functions 

Y J ord x J{x) = 0. (2.11) 

This result can be interpreted as some sort of conservation law for Coulombic- 
like "charges". This analogy is not superficial. It will be discussed in a separate 
publication. 

Eq.(2.11) is valid for both number and function fields and is known number- 
theoretically as the product formula 

n p H P =i- (2.i2) 

Taking log of both sides of Eq.(2.12) and using definitions of valuations one 
reobtains Eq.(2.11) back as required. From the discussion we had so far it is 
clear that the order and the degree have almost the same meaning. In particular, 

every element f(T) e F[T] has the form f(T)=a T n + ct x T"" 1 H V ct n and 

if «o 7^ it can be said that deg(/) = n. If / and g are non-zero polynomials, 
we have [42,43] 

deg(/.g) = dcg(/) + deg( 5 ) 

and 

dcg(/ + g) < max(dcg(/), dcg(g)). 

Thus, in view of the definitions given above, for the function fields the valuations 
are always non-archimedean. 

We would like to take advantage now of the fact that valuations can be 
introduced in many ways. In particular, the finite field F p with p elements is 
obtained as a quotient Z/pZ with p being a prime. If the coefficients a» of a 
polynomial f(T) of degree n belong to F p , then the quotient F p [T]//(T) = 
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GF(p") = F q forms the Galois field of q = p n elements (congruence classes), 
e.g. see Ref.[49], page 422, and Appendix A of this work. The advantage of 
working with finite fields lies in the fact that for each n there is an irreducible 
polynomial in F p [T], e.g. see Ref. [49], page 472, and Appendix A of this work. 
This means that if we use the irreducible polynomials only, then the degree can 
be used instead of order. To facilitate reader's understanding of these facts, 
a summary of relevant number theoretic results is provided in the Appendix 
A. In view of these results instead of Eq.(2.6) it is more advantageous to use 
|/| = q d °s(/) as valuation and it can be checked [12,42,47] that such defined 
valuation obeys axioms a), b) and Eq.(2.3) and, hence, is non-archimcdean. 



2.2 Zeta functions 

Definition 2. 3. The zeta function ( F (s) of the function field F q [T] is defined 
by 

Cf„ ( s ) — ^ TTp" i where / is the monic polynomial. (2-13) 

f£F q [T] '■'I 

Recall [46] that the monic polynomial f(T) is such for which a>o = 1. Using 
Appendix A we find that there are exactly q d monic polynomials of degree d in 
F,[T]. Hence, one obtains [43] 

deg(/)<d U 1 H 1 H 

and, accordingly, 

In view of Eq.(1.23) we obtain as well Cf„( s ) = Cq( s ~ !)• Next, we define zeta 
function over the rational function field K= F(T). For this, we need a notion 
of the projecive line P 1 (F) for the function fields. Since the usual projective 
line is just R U{oo}, evidently for the field F q it is K= F q U {oo}. In our case 
we have two types of polynomials, e.g. see Eq.s(2.8) and (2.10). Clearly, Cf, ( s ) 
is applicable for the first type. If we include the "prime at infinity", we obtain 

Ck(*) = ; \ \ 1 — • ( 2 - 16 ) 

Let now h(f) be some function (to be discussed explicitly later) defined on the 
set of monic polynomials. With help of this function we obtain the following 
Definition 2.4. The Dirichlet function Dh(s) associated with h(f) is given 

by 



oo 



D h (s)= ^ »>. = ^H(n)u n (2.17) 

/ monic n—0 
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where u = q s and 

H(n)= Y, M/)- (2-18) 

/ monic 
dog(/)<n 

The function H(n) should possess the multiplicative property, that is H(m)H(n) -■ 
H(mn). The Dirichlet function Dh,(s) sometimes is also called as zeta function 
Zk(u) associated with the field K [43]. One of the most fundamental achieve- 
ments of number theory is summarized in the following theorem 

Theorem 2.5. Let K be a function field in one variable T with a finite 
constant field F with q elements. Suppose that the genus of K is g. Then there 
is a polynomial Lk{u) £ Z[u] of degree 2g such that 

Zk(u) = - £ *ff (2.19) 
(1 - qu)(l - u) 

where L K (u) = (1 — c^u) with aia g+i = q, 1 < i < g. 

i=l 

The proof of this theorem, attributed to A.Weil, can be found in many places, 
e.g. see Refs[17,43-45]. 

For the sake of space we refer our readers to the literature for details of 
relation between the genus g and the field K. In short, we can embed f(T) € 
F q [T] into n-dimensional projective space P"(F 9 ) that is to convert it into the 
homogenous polynomial and to equate such polynomial to zero. The genus 
is an invariant of this homogenized polynomial with respect to some known 
transformations. All this is beautifully explained by Raoul Bott in Rcf.[50]. 

Combining Eq.s(2.17)-(2.19) produces 

oo 

Z K (u) = J2 H ( n ) un - ( 2 - 20 ) 

We shall need yet another interpretation of the same results. To facilitate mat- 
ters, in view of Eq.(2.15), following Rosen, Ref. [43], we write formally 

1 - qu 11 

where a; is number of monic irreducible polynomials of degree I . To obtain a; 
we need to take the logarithmic derivatives with respect to u of both sides of 
Eq.(2.21). Then, multiplying the result by u produces 

<l u = la i ul (2 22) 

1 - qu ^ 1 - u 1 ' [ ' ' 

H i=i 

Finally, upon expansion comparing the coefficients of u l produces the following 
result 

X>J=9" (2-23) 

l\n 
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which provides the missing link between a; and q. Looking at Eqs(2.15),(2.17),(2.20) 
and (2.21) one can write as well for general case of genus g similar result: 



Zk(u) = H(1-u 1 ) 



(2.24) 



1=1 



This time, however, the connection between a/ and q is more complicated. To 
obtain this connection we have to take the logarithm of both sides of Eq.(2.24). 
By expanding this logarithm in a power series in u we obtain, 



Here, by definition, N m is the number of zeros of f(T) in P"(F«j). For the 
projective curve of genus g associated with such polynomial, it can be shown 
[17,43-45] that 



For g = substitution of this result into Eq.(2.25) reproduces back Eq.(2.16) 
as required. It should be noted, however, that the results just presented arc 
only valid for the so called non-singular curves, that is for the Ricmann surfaces 
without nodes and holes. For singular surfaces (curves) the above formalism 
should be amended as discussed (albeit incompletely) in Ref. [51]. To bring all 
these results closer to those obtained in section 3 of our earlier work, Ref. [4] , it is 
helpful to provide the alternative (dynamical) interpretation of just obtained 
results. 

To this purpose, we begin with simple observations. Consider a map / : x — * 
y given by y = f(x). The number of the fixed points of this map can be easily 
obtained graphically by plotting together at the same plot the straight line y = x 
and y = /(x). Following Ref. [52] we define the global Lefshetz number £(/) as the 
intersection number /(A, graph(f)) where A is diagonal {x,x) in the direct 
product (x, y) G X x Y. Next, we need a notion of the local Lefshetz number. It 
can be easily understood if we recall how the linearization of dynamical systems 
is done. For instance, for two dimensional case we have f(x) = Ax + higher 
order terms where the dynamical matrix A can always be brought to the diagonal 
form and, hence, given by 



so that the local Lefshetz index Co(f) is given by sign [(ai — l)(c*2 — 1)] . For 
a source or a sink both eigenvalues have the same sign while for a saddle they 
have the opposite sign. Clearly, the local Lefshetz number is equal to the index 
of the vector field, e.g. see our earlier work, Ref. [6], for illustrations and further 
details. Hence, 




m— 1 



(2.25) 




(2.26) 





(2.27) 



f(x)=x 
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and this is yet another way of writing the Poincare-Hopf index theorem exten- 
sively used in our work, Ref.[6]. In this reference the index I(x) = C x (f) of 
singularity at point x was defined as 

I(x) = ±-lmf^ (2.28) 

c 

with C being a unit circle. This definition is valid for two dimensional case only. 
It is easy to extend this result to higher dimensions. Following Ref.[52] the local 
Lefshetz number can be identified as well with the degree (deg) of mapping <p : 

~ • f{z) - Z (2.29) 



More accurately, suppose z is an isolated fixed point of f(x) in R fc (any manifold 
locally is R fc ). If B is a small closed ball centered at z that contains no other 
fixed points, then the map ip is a smooth map defined by ip : dB — > S k ~ 1 . For 
k = 2 this is just mapping of one circle onto another (or to itself). The number 
of times one circle winds around another (the winding number at the point z) is 
deg z <p. In view of these remarks, such definition of the degree clearly coincides 
with Eq.(2.28) in two dimensions. Hence, I(z) = C z (f) = deg z (p. 

These results can be easily generalized now as follows. Following Milnor[53], 
let us consider sequence of iterates ,f 2 (z) = f(f(z)), f 3 (x) = f (f (f (z))) , etc. 
and their fixed points. Clearly, we will obtain the associated Lefshetz numbers 
£(/ 2 ), £(.f 3 ), and, accordingly, we can construct the zeta function 

lnZ f (u) = ^^u m (2.30) 

m— 1 

whose meaning should be clear. Obtained results are to be used in the next 
section. 



3 Variations on a theme by Fermat 

3.1 Circle maps, function fields and the Veneziano ampli- 
tude 

We would like to extend connections between the number theory and the theory 
of dynamical systems in order to reexamine results obtained earlier in Ref.[4]. 
Let us begin with the observation that every polynomial is completely defined by 
its coefficients a (Appendix A). Suppose now that these coefficients belong to the 
field F p and consider multiplication of such polynomials. This multiplication 
can be performed as follows. First, we remove the requirement that p is the 
prime number. Next, we introduce the space O^v (where N > 2 replaces p) 

Q N :={uj = (c , ci,C2, ...) | d € [0, 1, ...,N - 1] for i e Z} (3.1) 
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of one sided sequences of N symbols which is in one-to-one correspondence 
with the set of polynomials. Multiplication of polynomials simply replaces to 
by u'=(c , c x c' 2 , ...) with coefficients c- <G [0, 1,...,N — 1]. That is polynomial 
multiplication corresponds to some permutation in the sequence of c^s i.e. to the 
automorphism : Ojv — > ^jv- Such procedure can be computerized and serves as 
coding/decoding system (cryptography) for transmission of information [44,45]. 
A particular realization of such automorphism known as one sided TV-shift can 
be defined as follows 

CT#:ft£->fi£ ,o%w = <J = K,44...). (3.2) 

The upper script R stands for the " right" . Surely, one can consider the iterates 
of such defined shift, and, accordingly, one can look for the fixed (periodic) 
points (orbits) for such iterates. 

Definition 3.1. The periodic orbits for the one sided shift are the periodic 
sequences (cjv) w = w if and only if c m+n = c n . Every periodic sequence is 
uniquely determined by its basis Co, ci,C2, Cn-i ■ 

Using results of Appendix A we conclude that there are N n periodic se- 
quences. In particular, if N is the prime number we get p n periodic sequences. 
This is the number of irreducible monic polynomials as discussed in Section 2. 
Hence, the dynamical zeta function for the one sided shifts coincides with that 
given by Eq.(2.15). This result is in accord with that obtained in Ref.[54] ( e.g. 
see page 107, Eq.(3.1.5)) by slightly different set of arguments. Being armed 
with such result we are ready now to accomplish much more. 

To this purpose we would like to remind our readers some results from the 
Nielsen-Thurston theory of surface homeomorphisms. As we have discussed 
in our previous works, Refs.[4-6], 2+1 gravity can be considered essentially 
as physically reformulated Nielsen-Thurston theory. For convenience of our 
readers the key points of this theory are summarized in the Appendix B. In 
this subsection we will develop our formalism without interruption referring 
to the Appendix B whenever it is necessary. In particular, using Proposition 
B.2., it should be clear that any surface homeomorphism can be lifted to the 
unit disk model of H 2 . Then, the surface dynamics becomes dynamics of maps 
of the circle. Using Appendix B and Rcf.[55] it follows that it is sufficient 
to restrict ourself by the periodic maps. If this is the case, then the circle 
S^can be identified with F m = Z/mZ and, if we are interested in the invcrtiblc 
transformations, then the number m should be some prime, i.e. m = p. Hence, 
we arrive once again at the cyclotomic fixed point equation, Eq.(A3), written 
now as 

f:x p = x. (3.3) 
If we ignore the trivial solution x = 0, it can be rewritten as 

x v ~ x = 1 (3.4) 

and will possess p — 1 nontrivial fixed points, n times iteration of Eq.(3.3) 
produces 

/" : x p " = x (3.5) 
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and, again, it will have p n — 1 non trivial fixed points so that the Lcfshctz 
number associated with the degree of such mapping is C x {f n ) = p n — 1 in view 
of Eq.(3.4). This result coincides with that obtained by different methods in 
Ref.[54], Proposition 8.2.4. Use of this number in Eq.(2.30) produces at once 

( oo n _ i \ 

= exp{- ln(l - pi) + ln(l - t)} 
1-t 
1-pt' 

Finally, by replacing t with p~ s in Eq.(3.6) we obtain the p-adic version of the 
Veneziano amplitude, Eq.(1.16b), discussed in Section 1. 



3.2 Prom Brieskorn-Pham to Fermat and back 

Previous result, as good as it is, by itself provides no clues about its extension 
to the multiparticle amplitudes. In this subsection we are going to make the 
first step towards this goal. The rest of the paper depends crucially on this first 
step. 

To begin, consider a set of polynomials Pb-p of the type 



B-P 



: /(z) = +*«!+... + & (3.7) 



defined in C™ +1 with oo, a n being an n + 1 tuple of positive integers greater 
than or equal to 2. The polynomial f(z) will be considered as the map /: C rl+1 
— ► C. The Brieskorn-Pham (B-P) variety Vb(/) is defined now as [25-27] 

V B - P (f) C n+1 | /(z) = 0}. (3.8) 

This variety has singularity only at one point: z = 0. The Fermat variety is a 
special case of the B-P variety for which ao = ai = • • • = a n = a. 

Remark 3.2. The Fermat variety is the Calabi-Yau variety when a=n + 1, 
Ref.[56]. 

Corollary 3.3. Because of this observation it is clear that if the Fermat 
variety has something to do with physics, then the methods used in mirror 
symmetry calculations [57] in physics can be applied directly to the Fermat 
varieties as well. 

Fortunately this is the case as will be demonstrated in this subsection. Our 
goals are broader than just to establish this correspondence since we want to cre- 
ate an additional links with number and knot theories which, to our knowledge, 
are absent in the existing treatments [57]. 

Consider an intersection K=Ve_p(/)n S 2n+1 (this intersection is actually 
a knot or link as it will become apparent from the discussed below) of the B-P 
variety with the boundary of 2n + 2 dimensional ball B and define the mapping 

*) - (3-9) 
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It is clear that </>(z) is again the circle map of the type discussed in the previous 
subsection. This time, it is more than just this map since, actually, the mapping 
<j) : S 2n+1 \K -> S 1 defines a smooth fiber bundle Ref.[58], Lemma 19.10, with 
base a circle S' 1 and a fiber the complement of a knot /link K in S 2n+1 . Earlier, 
in Ref. [4], Section 3, we discussed in detail the mapping torus construction of 
3-manifolds fibering over the circle. We would like to argue here that for the 
Bricskorn variety associated with the polynomial f(z) = Zq° +Zj' we can recover 
back results obtained earlier. 

To this purpose let us recall the mapping torus construction first. Consider 
an orientation preserving surface homeomorphism h : S — > S. In the case of 
punctured (holed) torus this homeomorphism should respect the presence of 
a hole. The circumference of the hole is our base space <S 1 (e.g. read Ref. [59]) 
which is our knot K, since the knot is just a circle embedded into <S 3 .The Seifcrt 
surface of the knot is a fiber and the 3-manifold is just a fiber bundle constructed 
in a following way. Begin with the product 5x0 (the initial state) and Sh x 2n 
(the final state) so that for each point x <E S we have (x,0) and (h(x),2n) 
respectively. The interval / = (0, 2ir) can be closed (to form a circle S 1 ) by 
identifying and 2ir causing the identification: 

i : (x,0) = (h(x),2n). (3.10) 

The mapping torus Th fiber bundle is just the quotient 

T h = ^. (3.11) 

It is 3-manifold which fibers over the circle and is complementary to the fibered 
knot (lying at the boundary of such (cusped) 3-manifold as explained in Section 
5 of Ref. [4]). 

In the present case let us consider the monodromy map h(z) given by 

h t (z ,...,z n ) = (exp{it/a }z , ...,exp{it/a n }z n ). (3.12) 

If the point ZQ,...,z n belongs to the variety Vg_p(/) then, the identification 
analogous to that given by Eq.(3.10) takes place by construction. And, hence, 
the rest of the fiber bundle related arguments also goes through. Therefore, 
indeed, the inverse map 4>~ (z) produces the fiber bundle known as Milnor 
fibration. It is n + 2 dimensional manifold embedded into S 2n+1 . Let n — 1, then 
we are dealing with the usual knots embedded in S 3, and their complementary 
3— manifolds. It is of interest for us to reobtain the Alexander polynomial Ak 
for such knots. This task can be accomplished if the monodromy matrix M 
is known as it is explained in detail in our earlier work, Ref. [4]. Then, the 
polynomial A K is given by 

A K (i) = det(iE-M) (3.13) 

with E being the unit matrix. Although we had discussed in Ref. [4] how the 
matrix M can be obtained, here we need yet another way of obtaining this 
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matrix. To this purpose following Pham [25] and further explained by Milnor 
[26], let us look at Eq.(3.9) and select such fiber for which 

This can be achieved, for example, by requiring /(z) = 1 or 

z^+z^+---+z^ = l. (3.15) 

Remark 3.4. For n = 1 and a = a\ = m this is just the standard Fermat 
equation written in the affinc form. 

Following Pham again, let us replace z?* by ti so that we get Ym=o^ = 1- 

If we impose an extra requirement : ti > 0, then such coordinates become the 

barycentric coordinates for simplicial complexes used in algebraic topology for 

homological calculations [60]. It is clear, that the only choice for z^ to become 

j_ 

complex while still to obey Eq.(3.15) is to replace ti by Zi — t"* exp(±^ I! ) . 

Next, looking at Eq.s(3.12) and (3.15) let us consider the homeomorphism 
ta : ^ , i.e. 

ta(z) = (exp{i27r/a }zo, exp{i27r/a„}z„). (3.16) 

This homeomorphism can be represented as a composition of homeomorphisms 
acting on the homology basis: one generator for one basis element. The basis 
elements are just circles, more exactly, the quotients lo € {Z/oqZ, ...,Z/a„Z}. 
The generators r a are associated with rotations, i.e. 

27TZZ/ 

r ai (u) = cxp(± )lo = Cuo (3.17a) 

(li 

with v being in the range 1 < v < aj — 1. Hence, is just an eigenvalue of 
the rotation (homology) operator. The reader should consult Pham [25] and 
Milnor [26] for more details. For the case of torus we have just two elements in 
the homology basis so that we have 

[r Bl ®r„ 2 ](wi,w 2 ) =£V(wi,w 2 ). (3.17b) 

With help of above information the Alexander polynomial can be written now 
according to Milnor [26] (and Bricskorn [27] ) as 

A K (*)= 1] (*- w i w 2), (3.18) 

where all Greek letters denote the corresponding roots of unity. 

To get a feeling of this result, let us consider the simplest example of the 
trefoil knot T. It was discussed in our earlier work, Ref.[4], and it is a typical 
example of torus knot [59,61]. More specifically, it is 2,3-type of torus knot. 
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Therefore, for such knot we have two " homological" cyclotomic equations: oj 2 = 
1 and uj 3 = l.The first produces just one relevant root lo\ = —1, while the second 



The result coincides with earlier obtained, Eq.(3.9) of Ref.[4], as required. Ear- 
lier, in Ref.[4], we argued, in accord with Ref.[59], that there are only two fibered 
knots for Seifert surfaces of genus one. The figure 8 knot is the only knot other 
than trefoil which has punctured torus as its Seifert surface. The Alexander 
polynomial for the figure 8 knot cannot be obtained from the polynomial given 
by Eq.(3.18) since the eigenvalues of the monodromy matrix always have mod- 
ulus equal to one by construction [62]. We had explained in Ref.[4] that at 
least one of the eigenvalues of the monodromy matrix should be strictly larger 
than one in order to reproduce the figure 8 Alexander polynomial. The ques- 
tion arises: what happens if we use numbers other than 2 and 3 describing the 
trefoil ? This issue was also addressed by Milnor [26] and Bricskorn [27]. The 
resolution of the apparent paradox lies in considering the Alexander polynomial 
for links instead of single knots. But this should be done with some care since 
links may require multicomponent Alexander polynomial for their description. 
This topic is discussed further in section 5. In principle, the beauty of this 
approach to knots and links lies in the fact that it is not restricted to knots 
and links embedded in S 3 . Multidimensional analogs of knots and links can be 
readily considered, e.g. with help of the generalized Alexander polynomial 



where k = — n and each ik runs between < ik < a^.This more general case 
is relevant for calculations associated with multiparticle Veneziano amplitudes. 
It is discussed in Section 5. 

Being armed with these results we would like to rcobtain them now in another 
way. This will bring us new physical interpretation of mathematically known 
results. Let us begin with the Fermat equation written in its standard form 



For N = 2 this is just the Pythagorean equation. The problem of finding 
all integer solutions for this equation is known as problem about finding the 
Pythagorean triples. We had discussed this problem in Section 2 of Ref.[4]. 
This problem is completely solvable. For N > 2 the problem of proving that the 
above equation has only trivial solutions (the Fermat 's last theorem) was solved 
by Andrew Wiles only in 1995 [63]. Some physically relevant results dependent 
on his proof are discussed in Section 4.3. At the same time, for finite fields 
the above equation has many solutions. To realize that such Fermat equation 




A(t) = J | (t - w wi • • • w„) 



(3.20) 




(3.21) 
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has actually many solutions it is sufficient to look at solutions of cyclotomic 
equations (using results on Appendix A) 

(3.22a) 

With their help Eq.(3.21) can be replaced by the linear equation in the projective 
space CP 2 

x + y = z. (3.22b) 

Eq. (3.22b) represents a hyperplane in such space. This hyperplane can be em- 
bedded into the Grassmanian associated with such projective space. This is 
done with help of Pliicker embedding so that each hyperplane (that is the par- 
ticular solution of the Fermat's equation) produces a point in the Grassmanian. 
The situation in the present case becomes analogous to that we had discussed 
earlier in connection with the Witten-Kontsevich model, Ref.[64]. With each 
point in the Grassmanian associated particular Veneziano amplitude as we shall 
demonstrate below and in Section 5. Mathematically, such an amplitude is in- 
terpreted as one of the periods of the Fermat's curve. Section 5 provides 
generalization of this result to Fermat hypersurfaces. 

In his fundamental work, Ref. [65] , Andre Weil had made estimates of number 
of solutions of Eq.(3.21) in cyclotomic fields and came up with his famous 
zeta function, Eq.(2.19). Below, in the next subsection, we shall demonstrate 
additional physical uses of such zeta function. For the time being, we would 
like to focus attention of our readers on issues of immediate relevance. 

In particular, Eq.(3.21) is written in the so called projective form. In such 
form it represents the Fermat curve in the complex projective space PC 2 . It 
may be also convenient to write it in the affine form, i.e. 

x N + y N = l, (3.23) 

by dividing both sides of Eq.(3.21) by z N . For N > 2, following Gross [34], 
we denote such affine Fermat curve as F(N). Clearly, F(N) is just a special 
case of Eq.(3.15). This time, however, our treatment of Eq.(3.23) is going to be 
different. From the algebraic geometry it is known [66] that genus g of F(N) is 
equal to ^(N — 1)(N — 2). That is the Fermat curve is the Riemann surface of 
genus g . We associate with such surface its Jacobian J(N) and subsequently we 
compute, following Gross and Rohrlich [34], the periods of differential 1-forms 
(to be constructed momentarily) on J(N). Let us call a double of integers (r, s) 
admissible if 1 < r, s < N—l (e.g. compare with Eq.(3.17a)). To any admissible 
double we associate the differential 

(3.24) 

of the second kind "living" on F(N). If £ is a primitive N th root of unity, let 
A and B be the automorphisms of the curve given by the formulas 

A(x, y) = ((x, y), B(x, y)={x, Qy). (3.25) 
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x, y 



N 



z. 



Vr,, 



dx 

,N-1 



The group generated by A and B acts naturally on -q r s . The forms r\ r s are the 
eigenforms for this action, with eigenvalues £ r J+ sfe . That is 9 

(A^B k ) VrtS = C j+S %,s- (3-26) 

This result should be compared with Eq.(3.17b). Based on the result by Rohrlich 
(Appendix of Ref. [34]) it is possible to show that this similarity is not just 
coincidental. For the sake of space we only provide a sketch of an argument 
leaving details aside. These details are given in Refs[35, 67]. 

The affine version of the Format curve, Eq.(3.23), can be parametrized as 
follows 

x = t, y = (l-t N )N. (3.27) 

Such parametrization naturally limits the range of t to the segment [0, 1] . Hence, 
we have effectively the mapping -f : [0, 1] — > F(N) of the segment [0, 1] to the 
Fermat curve. By combining this result with Eq.(3.24) we obtain, 

i - /V-v-i_^_ = lf-H\-t N y^r d l = 1 B (- -) 

,%s ~i V V"- 1 J ( ' (1-tiV)^ 1 N B[ N>N>- 

(3.28) 

Let = a, jj = b, then we obtain the nonsymmetrized Veneziano amplitude, 
Eq.(1.2), 

_ r( )r(6) 

iVir ' s_ T(a + b)- 

Clearly, the symmetrized Veneziano amplitude can be obtained now without 
any problems. 

Remark 3. 5. The numbers a, b and c = 1 — a — b, can be made negative in 
accord with Veneziano formula, Eq. (1.6). The details are explained in Section 5. 

To obtain the symmetrized amplitude and also to make connection with 
knots and links ( in view of Eq.s(3.17b), (3.26)), we need to make few additional 
steps. Following Lang [35], let ( — exp(^p). Using the affine version of the 
Format curve, Eq.(3.23), it can be shown that the chain 

« = 7 - (1, C)7o + (C, C)7o - (C, l)7o. (3-30) 

where (1, Q, etc. have the same meaning as in Eq.(3.25) and k is a closed path, 
i.e. a cycle on the Fermat curve. When lifted to the Jacobian J(N) such a 
curve becomes one of the toral periods. Hence, using Eq.(3.30) and making 
trivial changes of variables one finds 



/ 



^ = (1-C r )(l-C S )^B(~). (3.31) 



Since 



9 See Section 5.3.2. for additional details 
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where t = N — r — s, we can rewrite Eq.(3.31) (up to a constant) in a manifestly 
symmetric form: 

/*,.,= -CK*r ((0 r, t ,r W [ sl „™»„^»„„]. (3.32) 

K 

Finally, by multiplying both sides by £2 W e obtain (up to constant again) 

= T(a)T(b)T(c) [sin 2a + sin 26 + sin 2c] . 

The r.h.s. of Eq.(3.33) looks almost the same as the Veneziano amplitude 
Eq.(1.12). It will have exactly the same particle spectrum. Surely, it can be 
made to look exactly the same if we use the symmetrized form of Eq.(3.28). 
Accordingly we obtain the following 

Corollary 3. 6. With accuracy up to a root of unity, the Veneziano ampli- 
tude is just a period of the Jacobian variety for the Fermat curve. 

The meaning of the phase factors ( r , ( s will be explained from another 
point of view in Section 5. 

Next, it can be shown [34,35] that 

V r ,s = J (A ] B k ) r? r;S = C 3+Sk J V r , s (3-34) 

and, therefore, the periods of the Fermat curve are eigenfunctions of the tensor 
product of the rotational (monodromy) operators, e.g. see Eq.s (3. 17a), (3. 17b), 
acting on the homology basis generated by closed n contours. From here we 
obtain yet another 

Corollary 3.7. Because the parameters r and s in Eq.(3.31) have the same 
meaning asp and q in Eq.(3.18) for the Alexander polynomial, the systematics of 
particle mass spectrum is in one-to one correspondence with the systematics of 
torus knot/links described by the Alexander polynomial, Eq.(3.18). More gen- 
eral polynomial, Eq.(3.20), describing multidimensional knots/links is discussed 
further in Section 5 10 . 

Relevance and connections of these results with knots/links will be discussed 
further in Section 5 devoted to calculation of the multiparticle Veneziano ampli- 
tudes. In the meantime, we would like to discuss some additional links between 
the number theory and physics by reinterpreting mathematical work of Koblitz 
and Gross, Ref. [68], in physical terms. 

3.3 Physical applications of Gross-Koblitz results 

In Section 2 we introduced Weil's zeta function, Eq.(2.19). From the associated 
discussion it should be clear that this function is effective measure of complexity 

10 In Section 5 we argue that, actually, it is more advantageous to relate the particle mass 
spectrum to the Hodge spectrum (defined in Section 5 and Appendix D) 
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of Riemann surfaces of genus g > 1 as compared to the "round sphere". This 
complexity is reflected in the polynomial L K {u) = Y\{1 — Q-iU) whose inverse 

i=l 

roots on carry all information about this complexity. The task of finding these 
roots was left unaccomplished in Section 2. Although the task of explicit finding 
of a^s can be rather complicated for arbitrary algebraic curve, for the Fermat 
curves this task can be brought to completion as was demonstrated by Gross 
and Koblitz, Ref.[68], see also [69]. In order to connect their results with the 
results discussed earlier in this paper let us consider the simplest nontrivial case 
of the elliptic curve, i.e. the case when g =1. Then Eq.(2.19) is reduced to 

7 , , {l-au){\ -qu/a) 

Zk(u) = — — — . (3.35 

(l-qu)(l-u) 

It should be clear that a and q/a are associated with homology basis for torus. 
They might me thought as two periods of the elliptic curve, more exactly two 
q-adic periods (as explained below in Section 4.3). To bring all this "down to 
earth" consider, following Landau and Lifshitz, Ref.[70], a simple problem about 
dynamics of heavy pendulum. The period T of such pendulum is determined 
essentially by the complete elliptic integral of the first kind : T = AK{k) 1 where 
k is known parameter and 

"2 

K(k) = [ dx (3.36) 

J v 1 — k sin x 
o 

Let, for example, k = l/y/2. Then, it can be shown [72] that 

to be compared with Eq.(3.28). 

Two questions immediately arise: a) is this a pure coincidence that mathe- 
matical expression for the period of the pendulum coincides with that obtained 
earlier for the Fermat curve? b) can heavy pendulum possess the second period 
and if it can what it means physically? To answer the first question we invoke 
the important result of Gerhard Frey [71] who discovered that the arithmetical 
properties of the elliptic curve Ea,b,c defined by the Weierstrass equation 

y 2 = x(x-A)(x + B) (3.38) 

are related to the arithmetic (diophantine) properties of the Fermat equation 

a p + W + c p = (3.39) 

which is reduced to equation 

A + B + C = (3.40) 
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if A = a p , B = 6 p and C = cP . This observation paved the way for proving the 
Fermat last theorem by Wiles [63]. To answer the second question(s) we follow 
Mark Kac (as described in Ref.[72], page 77). The argument goes as follows. 
Suppose we change time in Newton's equation for pendulum : t — > \J — It. Such 
substitution is equivalent to reversing of direction of the gravitational force and 
leads to a complementary periodic motion of the pendulum. That is the solution 
of the pendulum problem has not only its obvious real temporary period but 
also purely imaginary period as well. It is given explicitly by \^lK(k') where 
k' = \f\ — k 2 . For k = 1/ \/2 we get k — k' so that the ratio of two periods is just 
But is the Gaussian prime (e.g. read our earlier work, Ref.[4]) that 

is it plays the same role as q in Eq.(3.35) and, therefore, indeed, our conjecture 
about the meaning of a's in Eq.(3.35) is apparently correct. Gross and Koblitz 
went much further in their analysis. We believe, that it is worthwhile to discuss 
their results in some detail since there is some physics associated with them. 

To accomplish this task we need to introduce some definitions. For more 
details reader is urged to read pedagogically written book by Winnie Li [17]. 
Let F g be a finite field. For a eF 9 define the additive character ^£> a of the 
field F q as follows: ip a (x) — ip(ax) \/x <GF 9 .The multiplicative character \{ x ) 
is determined by the cyclotomic equation, Eq.(A.5). It is just one of the q- 
th nontrivial roots of unity and there are q — 1 of them. With help of these 
characters one can construct the Gaussian and the Jacobi sums. In particular, 
the Gauss sum is defined as 

9(X>1>) = - X(x)ip(x) (3.41) 

where F* =F g \ 0, while the Jacobi sum is defined as 

J(Xi,Xt) = - E Xi(*)x 2 (1-z)- (3-42) 

The Gauss sum is analog in F g of the gamma function 

oo 

/dr 
x°e-*- (3.43) 

o 

since this is "sum" (i.e. integral with respect to its Haar measure — ) over the 
multiplicative character x s and the additive character e~ x . The Jacobi sum is 
analog of Euler's beta function since 

J(Xi,X 2 ) = g{X ^ )g{X ^ ] . (3.44) 

g{xiX2,i>) 

The p-adic gamma function (not to be confused with Eq.(1.21)) T p (n) for an 
integer n > 2 is defined by 

n-l 

r p (n) = (-l)" I] 3- (3-45) 

(p,j)=i 
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It possesses almost standard properties, e.g.r p (x + l) = j p^?^f ^ X ^ ~ p ■ 

Koblitz had demonstrated in Ref.[69] that for the Fermat curve, Eq.(3.23), the 
zeta function of Weil is given by [see also [111], ch-r 6], 

Z(F(n):u\F q ) = — — L— — JJ (l-J(x' >X » (3.46) 

V q A ' l<r,s<JV-l 



In addition, he shows that 

Jt\'.\n= ~'l N ;z\ N ' (3-47) 



■^ p(iy)^p(jy) 
r c r +- s ^ 

P v JV / 



to be compared with Eq.(3.28). From here it follows that Eq.(3.47) is the p-adic 
analogue of the 4 particle Veneziano amplitude. In proving the Theorem 4.13 
of Rcf . [68] Gross and Koblitz show that knowledge of the p-adic beta function 
allows to restore the full beta function, i.e. the full Veneziano amplitude and 
vice versa. This observation provides the strongest support to the conjecture of 
Raoul Bott who said, Ref.[50], page 159, that "I would not be too surprised if 
discrete mod p mathematics and the p-adic numbers would eventually be of use 
in the building of models for very small phenomena" . 

Remark 3. 8. It is interesting to notice at this point that if instead of the 
complete elliptic integral, Eq.(3.36), we would consider an indefinite elliptic 
integral, e.g. 

dP 



I 



47> 3 - g 2 V - 53 

then, its inverse is the stationary solution of the KdV equation, i.e. 

(V) 2 = 4T 3 - g 2 T - g 3 . 

Mulase [72] had shown how one can easily obtain the time-dependent solution 
from the time-independent one. Hence, the dual of KdV leading to the Virasoro 
algebra (as we had discussed in our earlier work, Ref.[73]) can be obtained as 
well. 

Plausible as they are, these results can be considerably extended strength- 
ened and generalized beyond the existing string theory formalism. Steps in these 
directions are presented in the following sections. 

4 The Chowla-Selberg formula, the Fermat 's last 
theorem, the conformal field theory and the 
Veneziano amplitude 

4.1 Statement of the problem 
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Originally, the Chowla-Selberg formula was announced in 1949, Rcf. [74] , but was 
left practically unnoticed. This caused the authors to provide more details. The 
revised and extended paper received the widespread attention among mathe- 
maticians but was published only in 1967 [31]. Ironically, the Veneziano's paper 
had appeared in 1968 [1] while the Annals of Mathematics paper by Ramachan- 
dra [75] published in 1964 containing detailed derivation of the Chowla-Selberg 
main results (to appear only in 1967 !) also was left largely unnoticed. 

In this work we would like to make connections between these historical 
papers. We begin with the description of the problem Chowla and Selberg 
wanted to solve. Looking at Eq.s(3.36) and (3.37) they wanted to know if the 
result, Eq.(3.37), is mere coincidence or if it is general rule. To make things more 
interesting, let us recall that the classical theory of elliptic functions provides 
us with the following results. Periods K and K' are given through the elliptic 
theta functions [76] 

K=\6l{r), (4.1) 
K' = \el{- 1 -), (4.2) 



while r is determined implicitly through 

,2 _ et(r) 



9i(r) 

with 



(4.3) 



9 2 = 2q*H Q Hf and 63 = H H% (4.4) 
provided that q = exp(7rir) and functions H , Hi and H 2 are given by 

00 00 00 

H =i[(i-q 2k ) ,ff 1 =na+9 2fe ) ,ff 2 =n( i +^ 1 )- ( 4 - 5 ) 

fe=l k=l k=l 

Based on these results, it seems like a hopeless task to calculate either of periods 
exactly. Needless to say the connection between Eqs.(3.36) and (3.37) looks 
totally mysterious. Nevertheless, Chowla and Selberg had demonstrated that 
this connection is not an accident. Under certain conditions to be discussed 
below, the representation of periods as products/ratios of gamma functions is 
the only possibility. In view of the Veneziano formula, Eq.(1.2) (or (1.6)), the 
conditions under which such representation is possible certainly is of physical 
interest. 

Traditionally, calculation of the periods is done with help of the hyper- 
geometric function as we had discussed earlier in Ref.[77]. It is appropriate 
to remind some aspects of this connection now. We have K=^F(^, |, 1; fc 2 ), 
where the hypergeometric function is obtained as solution of the hypergeometric 
equation 

z((l - z)F" + [c - (a + b + l)z]F' - abF = 0. (4.6) 
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Since both K and K' can be obtained as solutions of Eq.(4.6) we can form their 
ratio 

m = Mt) = kM> (4J) 

where yiand t)2 are solutions of the Fuchsian -type equation 

y" + \{w,T}y = (4.8) 
with the Schwarzian derivative {w,t} known to be as 

(4.9) 

with w = w(t),w' = 4^, etc. Clearly, everybody who is familiar with string 
and conformal field theories will recognize at this moment connections between 
mathematics and physics. We are not going to develop such connections in 
this paper nevertheless. The only reason we have brought the hypcrgcometric 
equation to the attention of our readers is to emphasize that this is equation 
for the periods of the elliptic curve. As it is well known, the Knizhnik- 
Zamolodchikov (K-Z) equations are basic for all conformal statistical mechanical 
models [78-80] and are effectively reducible to the equations of hypergeometric 
type. Hence, all results of conformal field theories will remain unchanged in 
what follows 11 . In particular, the expression for the nearest neighbor spin-spin 
correlator for the Ising model involves -F(|, 5, 1; k 2 ) [81], page 69 . But, much 
more can be actually accomplished thanks to results obtained by Chowla and 
Selberg. 

To begin, we would like to remind our readers about the partition function 
for the free bosons on the torus described, for example, in Ref.[82], pages 340- 
344. The case of free fermions (Ising model) technically is almost the same. 
Indeed, in the first place the whole computation depends upon calculation of 
the sum 

G{s) = V , - , 2 „ (4.10) 



m,n |m + nr| 



while in the second, of the sum 



m.n \m + nr+ (h + vt)\ 

The purpose of the entire calculation in both cases lies in explicitly obtaining 
G"(0) and G'^ „(0) where the prime means differentiation with respect to s vari- 
able. It is clear that in both cases obtained expressions are r— dependent. They 
also should be invariant with respect to modular transformations: r — > r + 1 

11 Strong additional support of this claim is presented in Section 5.3.2 below. 
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and t — ► — - . It is well known, e.g. see our earlier work, Rcf.[4], that all 

T 

transformations of the type 

r> = ^± (4.12) 

with a, b, c, d being some integers subject to restriction ad — be = 1 can be ob- 
tained by successive applications of the above two. In the same work we also 
explained what these transformations mean topologically. In view of generaliza- 
tions which follow we need to present these results from the different perspective 
now. 



4.2 Complex multiplication 



To this purpose it is helpful to recall that, according to the theorem by Jacobi 
[76], the modular lattice L, 

L = Tiuoi + Zw 2 or, symbolically, L=[oj\,uj2\, (4-13) 

has periods such that r = Im — > 0,that is one of the periods should be 

complex. Eq.(4.13) looks the same as Eq.(A.2) and, actually, has the same 
meaning. Moreover, from the algebraic point of view, L is just an ideal. In 
the Appendix B of our work, Rcf.[4], we explained the notion of an ideal. For 
the sake of uninterrupted reading we would like to recall few things from that 
reference. 

Suppose we have some set / and another set L,then Va, /3 £ I and V£ G L, 
we have 

a + (3 G I (module property), a£ G L (ideal property). (4-14) 

In the present case, we have a and (3 € Z and L=[wi,U2]- Clearly, if u)\ and 
u>2 constitute the basic parallelogram of such lattice then, any point x=awi + 
(3iV2 also belongs to the lattice. Surely, for a given lattice we can make many 
sublattices using equivalence relation x = y (mod m) , that is x — y = mk for 
some integer k. Specifically, if we initially had lattice L=[u)\, U2] we can consider 
lattice £'=[u/ 1; u' 2 ] where, in view of the results of the Appendix A, we have to 
demand 

u> 2 = auj 2 + bu>i (4-15) 
oil = ojj 2 + dui. 

Eq.(4.12) is obtainable from these two equations by forming their ratio. We 
refrain from doing this for the moment since we want to consider the possibility 
u[ = au) 1 ,u> 2 = otuJ 2 - That is we want to study the eigenvalue problem for the 
matrix A given by 




(4.16) 
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According to Eq.(B.3) of Appendix B the problem of finding the eigenvalues 
of matrix A is reduced to the problem of finding eigenvalues of the quadratic 
equation 

a 2 -trA a + detA = 0. (4.17) 

This equation should be supplemented with some physical conditions in order 
to make sense of its solutions. In particular, let a(L) denote the area of the 
period parallelogram associated with the lattice L. Evidently, [83], 

a(L) = ^ |wiw 2 - w 2 wi| , (4-18) 

where the overbar means the complex conjugation. In particular, for the lattice 
L T = [1, r] Eq.(4.18) produces a{L) = |Im(r)| . Now, if we rescale oj's, we obtain 
a(aL) — \a\ a{L). If we require that, upon such rescaling, the area a(L) remains 
unchanged this leaves us with the only one option: \a\ =1. This result should 
be dealt with in connection with Eq.(4.17). This leaves us with not too many 
choices. One option is to have trA = thus producing rather trivial result: 
a = ±1. The rest of options comes from the standard number-theoretic results 
about units in the quadratic fields as discussed for example by Hardy and Wright 
[84] (see also Appendix B of our earlier work, Ref.[4]). Hence, in addition to ±1, 
we obtain also a = ±i and a = ±p, a = ±p 2 with p = |(— 1 + >/— 3). The last 
result we had encountered earlier in connection with the Alexander polynomial 
for trefoil knot, Ref.[4], Section 3. These are the only options available for the 
quadratic fields. In the light of these results the following definition is useful: 

Definition 4.1. Two complex tori T = C/L and T" = C/L' are isomorphic 
(respectively isogenous ) provided that there is an isomorphism (respectively 
isogeny) a : T — > T' induced by multiplication a : C — > C for a e C with 
aL = L' (respectively aL C L'). 

From this definition it follows that isogeny is an equivalence relation because 
L'is homothetic to L. By establishing such relation we are effectively embed- 
ding our tori into the projective space. Homothetic lattices are associated with 
isomorphic elliptic curves [85]. This feature, known in literature as complex 
multiplication ( CM) [85] , admits generalization to Riemann surfaces of higher 
genus [18]. 

It is important to realize that in general an arbitrary torus T may not possess 
CM. In the traditional (and, hence, more general) case there is no need to require 
uii = aoj 1 ,Lu' 2 = au 2 - In this case Eq.s(4.15) are reduced to Eq.(4.12), provided 
that ad— be = 1. Eq.(4.12) might be interpreted as equation representing motion 
in the Teichmiiller space of the punctured torus: Different initial points in the 
moduli space for T will produce different orbits in the Teichmiiller space of the 
torus. All this is explained in our earlier work, Ref. [4]. 

Let us now go back to Eq.s (4.15). The second of these equations can be 
rewritten as a = cr + d and, since r is a complex number, a should be a complex 
number as well. 

Remark 4.2. At this point we would like to notice that if r would be real, 
the multiplier a becomes real as well. Then, instead of CM theory one would 
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have real multiplication (RM) theory. Such theory would make no sense from 
the point of view of traditional algebraic geometry but would make sense from 
the point of view of non commutative geometry of Connes according to the 
recent papers by Manin [86,87]. In our work such possibility is not considered 
since CM is sufficient for description of physically relevant phenomena. Surely 
we do not exclude the possibility that RM might have some physical significance 
and leave this option open for further study. 

Going back to the subject of CM, we notice, that the first of Eq.s(4.15) can 
be written as (ct + d)r = ar + 6. This result can also be written in the form of 
Eq.(4.12), i.e. 

ar + b . A . 

r = —-j. (4.19) 

This is equation for the geodesic in H as we had discussed in Section 2 of Ref. [4] . 
Explicitly written, it looks like 

ct 2 + {d-a)T -6 = 0. (4.20) 

It can be equivalently rewritten as 

(cr) t = b + {d - o)t. (4.21) 

If we denote ct as (3 then, the above equation can be presented as (3uj2 = 
boji + (d — a)uj2 and, in addition, we have I3ui\ = clu 2 . These equations are of 
the same kind as Eq.s(4.15) as required. Moreover, using the same equations 
and eliminating t we obtain the following equation for a : 

a 2 - (a + d)a + ad - be = (4.22) 

in agreement with Eq.(4.17) or, taking into account that ad — be = 1, we obtain 

a 2 -{a + d)a + 1 = 0. (4.23) 

This is equation for an integer in the quadratic field (Appendix B, Ref. [4]). 
But earlier we have obtained a = ct + d and, hence, such integer must belong 
to the imaginary quadratic field. Therefore, r also belongs to the imaginary 
quadratic field but is not an integer in general since r = (a — d)/c. As we 
had discussed in Ref. [4], the imaginary quadratic fields are being character- 
ized by their class numbers h(>/—d). These numbers directly associated with the 
unique factorization property for the ideals of these fields. Every number x 
in imaginary quadratic field can be presented as h = a + b^Z—d, where a and 
6 e Q. Unique factorization producing class number h(\/—d) = 1 is possible if 
— d = 1,2,3,7,11,19,43,67 and 163. Sometimes, instead of 1 and 2 numbers 4 
and 8 being used but, surely, these produce the same results. We mention this 
fact for the sake of comparison with the existing literature. If we would like to 
restrict ourself by these numbers, the parameter r in Eq.s(4.10) and (4.11) can- 
not be arbitrary anymore. It should be discrete. Nevertheless, because 
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we can vary integers d and c in the equation r = (a — d)/ c,this discreteness does 
not lead to just one value for r. 

Now we would like to illustrate all these statements by simple but important 
examples. These examples serve to underscore the differences between what 
is known so far in the conformal field and string theories and new elements 
which, in our opinion, should be included. Inclusion of these new elements will 
enable us to use powerful number-theoretic methods in the conformal and string 
theories 12 . 

Going back to Eq.(4.23) we obtain the following expressions for a's : 



Oll,2 



a + d 1 
2 2 



(a + df - 4. 



(4.24) 



Since a should be an integer in the imaginary quadratic field we are left with 
the following set of options: 

a) a = d = 0, thus producing oi ; 2 = ±i; 

b) a = ±1, 6 = (or a = 0, b = ±1), thus producing ot\ t 2 = |(±1 ± V~ 3) 

c) a — d — ±l,thus producing 0:1,2 = ±1. 

We have obtained all these results earlier and, surely, we would not reproduce 
them once again should these be the only options. Fortunately, there are other 
less trivial options leading to far reaching consequences. 



Consider a(L') = 



where u>[ and u' 2 are given by the l.h.s. 

of Eq.(4.15). This expression can be rearranged with the help of Eq.(4.15) as 
follows: 



a(L') = 



1 

2 
1 

2 

n _ 
- l w i w 2 



■ W2 w 'l 



(4.25) 



\(ad — bc)(LO\uj2 



W 2 Wl 



- W2Wl)| 

na(L) 



provided that ad — be = n. We would like to argue now that such relation is 
legitimate. Indeed, the relation ad — be — 1 comes as result of the requirement 
for the inverse of the matrix A, Eq.(4.16), should have only integer coefficients 
[88]. This is a legitimate requirement as long as we are dealing with the same 
lattice and just trying to change the basis. But if we are dealing with sublattices, 
relation ad — be = 1 can be relaxed to ad— be = n. This can be rigorously proven, 
e.g. see Rcf. [89], ch-rl, and leads to a very deep results summarized in the 
Appendix C. The scaling result a(aL) = \a\ a(L) discussed earlier can now be 
rewritten as a(L') = a(ctL) = na(L). This suggests that if L is sublattice of L' ', 
then these sublattices are homothetic and therefore are isogenic and, hence, are 
equivalent. Accordingly, Eq.(4.22) should be amended leading to 



011, 2 



d 



2 ± l -^(a + dy-An. 



(4.26) 



12 Usefulness of these "new elements" to conformal and "string" theories can be seen from 
both: the references provided at the end of Section 1.3 and from actual developments presented 
in Sections 5.2. and 5.3. 
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Thus, indeed, the complex multiplication is associated with integers in the imag- 
inary quadratic field. From the above discussion it follows that there are only 9 
imaginary quadratic fields whose ring of integers has class number /i(v / — d) = 1. 
Actually, more can be said. Following Silverman, Ref.[85], we notice that the 
j(L) invariant of the elliptic curve (modular form of weight zero, e.g. see Eq.(5.8) 
below) is just some rational number e Q provided that j{L) characterizes the 
equivalence classes of elliptic curves with CM whose class number is 1. Moreover 
j(L) e Q only in 3 other cases so that there are exactly 13 classes of elliptic 
curves with complex multiplication. Such curves are known as Weil curves. 
They possess remarkable properties to be briefly discussed in the next subsec- 
tion. It should be noted, however, that we are concerned about h — 1 case only 
because it is technically much simpler as compared to other quadratic fields for 
which h > 1. Altogether, quadratic imaginary fields are much simpler to handle 
than the real quadratic fields [90]. 

4.3 Fermat's last theorem, periods of elliptic curves and 
Veneziano amplitudes 

Zcta function of the elliptic curve Zk{u) is given by Eq.(3.35). Now we would 
like to connect this function with modular forms of Appendix C. According 
to famous Taniyama-Shimura (T-S) conjecture (now proven, thanks to works 
of Wiles [63] and Taylor- Wiles [91]) every elliptic curve over Q is mod- 
ular. Naturally, we would like to explain what this notion means and how it 
is connected with Zk{u). Since the results to be discussed depend crucially on 
(now proven) T — S conjecture and since the proof of this conjecture effectively 
implies, according to works by Frey [71] and Ribet [92], the Fermat's last the- 
orem, e.g. see Eq.s(3.38)-(3.40), it happens, that the mathematics related 
to Veneziano amplitudes is also related to the Fermat's last theo- 
rem. Accordingly, the multiparticle generalization of these amplitudes, to be 
discussed in Section 5, is dependent to some extent also on the validity of the 
Fermat's last theorem. 

Let £ be an elliptic curve over Q. Analytically, it can be written in the so 
called (affine) Weierstrass form as 

y 2 + a\xy + a 3 y = x 3 + a 2 x 2 + a^x + a 6 (4.27) 

with coefficients in Q. Let A(ai, ae) be the discriminant of the elliptic curve 
£ [83,85] . It carries information about the smoothness of £. If A = the curve 
is singular, that is it may contain a node or a cusp. Suppose that our curve is 
nonsingular in the field Q. Suppose as well that we would like to know if such 
curve will still remain nonsingular if we reduce it, that is, if we would like to 
consider instead of Eq.(4.27) its F p analog (Appendix A) 

y 2 + a\xy + a 3 y = x 3 + a 2 x 2 + a^x + a 6 ( modp) (4.28) 

Remark 4. 3. Such reduction allows us to use without change all results of 
Section 2.2. 
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Surely, to do so we need to get rid of denominators in Eq.(4.27) so that all 
coefficients belong to Z. Upon reduction the curve £ p may become singular 
(even though it was not singular initially!). This happens when A(ai, a 6 ) = 
O(modp) and represents the case of bad reduction. Otherwise the reduction is 
called good. Clearly, £ has bad reduction only for a finite number of primes. 
The bad reduction can be measured by means of the conductor. The conductor 
is a number A such that 

Cond{£) = N = Y[p f ". (4.29) 

v 

Here f p = if p \ A and f p > 1 othervice. Let Ai be the number of solutions 
of Eq.(4.28) in the projective space (e.g. see Eq.(2.26)) then, the Hasse-Weil L 
function of £ can be defined as 

L{E, s) = Y[L p (p- s ), (4.30) 



where, in the case of good reduction, [L p (p s )] = (l—pu)(l—u)Z p (u) | u=p - s with 
Z p (u) defined by Eq.(3.35). In general case we have as well 

L p (p- S ) = 1 - a p p- s + ^y- 2s (4.31) 

with a p = 1 +p — Ni(p) and ip(p) = or 1 depending upon p dividing or not 
Cond{£). Hence, in general we can write 

1(g " )= S^g ,4 - 32) 

Now it is time to connect this result with the modular forms discussed in the 
Appendix C. More exactly, we need to make connections with the cusp forms 
of weight 2. The cusp form /(t) is defined by its q expansion: /(t) = 
S^Li c{n)q n ', q — exp(27rzr). Using the properties of Hecke operators (Appendix 
C) and, in particular, Eq.s(C.9) and (C.16) (for k=2), we obtain (upon acting 
on /(t)) the following set of recursions for c(n)' s : 

c(p e )c(p) = c(p e+1 ) + pcipf- 1 ) for p prime, p \ N, (4.33a) 

c(p e ) = [c(p)] e , for p prime and p \ A, and (4.33b) 
c(m)c(n) — c{mn). (4.33c) 
The Mellin transform of the cusp form can be defined now as 

oo ^ 

F(s) = /^ / ^)r-A( / , S ):^(2.)- s r( S )^^ (4.34) 

{ s 

EE (27T)- S r(.s)L(/, S ). 
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It can be easily demonstrated, e.g. see Ref. [93], that L(f, s) = L(£, s) provided 
that the cusp form (of weight k = 2) is also an eigenform for the involution 
operator Wn defined by 

W N f(r) = ±NT~ 2 f(j±). (4.35) 

More exactly, it is expected that the cusp form of weight 2 is an eigenform of 
two operators T 2 (n) and Wn ■ 

T 2 (n)f(T) = c(n)f(T), (4.36a) 

WnJ{t) = V f(r) (4.36b) 

with r/ = ±l(fermionic property). If /(r) lies in one of the two eigenspaces for 
Wn and satisfies Eq. (4.36a) then, L(f, s) = L(£, s). 

Definition 4.4. An elliptic curve over Q is said to be modular if there exist 
a cusp form of weight 2 on Sn (Appendix C) for some N such that L(f, s) = 
L(£,s). 

Remark 4.5. This identification means that a' p s in Eq.(4.32) should be 
replaced by c(p)'s 

Shimura and Taniyama conjectured that every elliptic curve over Q is modu- 
lar and Wiles [63] proved correctness this conjecture. Now it remains to explain 
what all these results have to do with the Veneziano amplitudes. We shall 
provide the answer based on results of the Eichler- Shimura theory (surely only 
those aspects that are relevant to our immediate needs) and keeping in mind 
further uses of the Chowla-Selberg formula (to be discussed later in Subsection 
4.5.). In addition, the entire Section 5 is devoted to elaborations of these results 
from yet another perspective. 

We begin with observation that the combination f(r)d,T is invariant with 
respect to action of modular transformations which belong to Sn- Define now 
a "path integral" 

7(to) 

<£/(7)= J WW. (4.37) 

TO 

This integral is actually independent of t for 7 e Sn- If 7 € Sn is elliptic 
(|ir7| < 2) or parabolic (\trj\ = 2) it can be shown [93] that $7(7) = 0. This 
is consistent with Eq.(4.20) which is expected to have two solutions, that is to 
produce the geodesic generated by the hyperbolic transformations. Hence, only 
hyperbolic transformations should be taken into account when calculating 
Since <f> f (7) is independent of To this is possible only if $ f (7) describes closed 
geodesic on the corresponding Riemann surface R = H/<Sn- But closed geodesic 
on the Riemann surface R corresponds to the period (e.g. see Eq.(3.31)) when 
it is lifted to the Jacobian J(R) of R. Surely, if R has more than one cusp form, 
there will be more than one period. Since these cusp forms are independent of 
each other, the corresponding periods are also independent. 
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Remark 4. 6. It is remarkable that such periods are defined on surfaces with 
cusps. In the case of dynamics of 2+1 gravity (equivalent to the dynamics of 
train tracks) the simplest non trivial Riemann surface on which the train track 
dynamics can be studied is represented by once punctured torus [4]. The sim- 
plest non trivial correlator for the Kontsevich-Witten model of 2 dimensional 
topological gravity also involves integration over the moduli space of once punc- 
tured torus [4]. 

The intersection pairing form < C , lo > can be defined now as 

f(£)d£ = J oj =< C,uj >= 0. (4.38) 
c 

This form can be easily defined for the Riemann surfaces of higher genus. Many 
interesting results related to Selberg trace formula, dynamics of hyperbolic flows, 
etc. can be found in Refs. [94,95] and references therein. In this work we are 
concerned with different issues however. In particular, we would like to explain 
why a and q/a defined in Eq.(3.35) can be called "periods" , more exactly, p-adic 
periods. We also would like to connect these p-adic periods with Q. 

Ignoring for the moment bad reduction terms in Eq.(4.32) it is useful to 
rewrite Eq.(3.35) in the following equivalent form 

i + iNp - p -i )p - +p L-*> = Cp(*)C P (*-i) , 39) 

(P ' (I-P-'XI-P 1 -*) - L p (s) (4 ' 39) 

where, according to Eq.(4.31), N p —p— 1 = — c(p) (see Remark 4.5). The global 
zeta function ((s) is obtained now as 

h = TT C P O0C P (s-i) = COOC(s-i) 
tU p L P (s) ~ L(£,8) 

with _ 

L(£, s) =]][1 + (N p -p- l)p- s + p 1 - 2 *] - 1 . (4.40) 
v 

Let s = 1 in Eq.(4.40). Then, we obtain a very remarkable result, Ref.[96], 

w,D=n(£). («i) 

Moreover, using Eq.(4.34), we obtain as well 

ioo 

L(£, 1) = -2ni J f{x)dx. (4.42) 
o 

Next, following Manin, Ref.[97], we would like to combine Eq.s(4.35), (4.36b) 
in order to get 

L(£, 1) = (1 ± 1) V ^ exp(-27rn/AT) (4.42) 
( n 

n=l 
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so that if c(n) are known, L{£ , 1) can be calculated numerically very efficiently. 
Next, the function L(f,s), Eq.(4.34), actually coincides with the Hecke zeta 
function L(s, \) with Grossencharakter x- It is known, that such function has 
the following presentation and the associated with it Euler product 

L ^ X) = E jM; = - x(p)N(p)-')-\ (4-43) 

a ^ ' p 

where the summation is taking place over all nonzero integral ideals with N(a) 
being the norm of the ideal (e.g. see Ref.[98] for a quick introduction to these 
concepts) while the product is taken over all prime ideals. In Rcfs.[99], ch-r.18, 
and [100], ch-r.8, it is shown (using some examples) that, actually, L(s,x) = 
L(f,s). Next, given this fact one needs to invoke the famous theorem by Weil 
[101] which identifies Hecke's Grossencharakters with the Jacobi sums, e.g. see 
Eq.(3.42). Once this is done, Eqs.(3.35),(3.38), (3.46), (3.47) and (4.31) should 
be used along with Definition 4.4. and Remark 4.5. in order to claim that, 

a(p) = a{r) + a(s) ~ J( X r , X s ) = ^ff^ = ^ • (4-44) 

Since according to Eq.(2.19) a ■ a = p,we obtain as well, a — p/a. Because 
the r.h.s. of Eq.(4.44) is a p-adic period, the l.h.s., that is a(r) (or a(p) ) is 
also a p-adic period. This proves the claim we made after Eq.(3.35) about the 
physical nature of a. Moreover, this proves the claim made earlier that to reach 
the conclusion about a the last Fermat's theorem should be used. 

Remark 4. 7. The connection between a's and the p-adic periods was first 
made by Dwork. Good and readable summary of this and related papers by 
Dwork is given in Rcf. [102]. Relevant results are also summarized briefly in the 
book by Koblitz , Ref . [69] . The arguments presented above differ however from 
those of Dwork. They rely on the fundamental result of Andre Weil [101]. 

We would like now to connect these facts with Eq.(4.41) and also with 
Eq.s(1.23), (1-24). To do this let us consider, following Birch and Swinnerton- 
Dycr, Rcf.[96], (and also Cassels, Ref. [103]) the product 

n^p- T (4-45) 

v 

where r denotes the Tamagawa number. Calculation of Tamagawa numbers is 
described for example in the book by Weil [104 ] 13 . Nevertheless, the problem 
of finding r in Eq.(4.45) is not completely solved to our knowledge. By analogy 
with Eq.(1.20), the above result remains unchanged if we replace oj p by |a| p Wp 
since n \a\ p = 1. That is the product, Eq.(4.45), is invariant with respect to 
p 

some scale changes. This is important in view of Eq.(4.44). The product in 

13 More recent useful reference is the monograph by Maclachlan and Reid, Ref. [105], where 
the Tamagawa numbers are discussed along with many other relevant facts. 
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Eq.(4.45) is taken over all places, including oo, therefore, 

r ~ Q Yl w p ( 4 - 46 ) 

finite p 

Surely, in view of Eq.(4.44) we expect f2 = B(r, s) with B(r, s) being the usual 
beta function that is the Veneziano amplitude. Eq(4.46) is the exact analog 
of Eq.(1.25) of the p-adic string theory [11]. To make all this connected with 
the number theory we have to look at Eq.(4.38) and to assume that 

J to «=* J iu p . (4.47) 

C Cp 

More accurately, this suggests that J £u p c=t c(p) and, since N p — p — 1 = —c(p) 

(according to Eq.s(4.31),(4.39)), one can divide both sides by p in order to 
obtain 

Np = p + 1 - c(p) 
p p 

In view of Eq.s(4.41),(4.44), we finally obtain 

-l 

(4.48) 

with the product being taken again over all places, including oo. The ramifi- 
cations of this identity lead to the so called Birch-Swinnerton-Dyer conjecture 
(still unproven). It can be stated formally as 

\im(s-l)- r L(£,s) =CCl TT oj p (4.49) 

finite p 

with C and r being some constants (details can be found in Ref. [85] ). Eval- 
uation of the constant C is the major stumbling block in proving the conjec- 
ture. The constant r is assumed to be known. Since this conjecture involves 
the Veneziano amplitude f2, Eq.(4.49) provides yet another link between the 
"string" and number theories. We would like to discuss now links between the 
number theory and the conformal field theories. 

4.4 Number-theoretic formulation of statistical mechani- 
cal models of conformal field theory 

The idea about deep links between the number theory and statistical mechanics 
is not new. It was discussed, for example, in the review paper by Tracy [106]. 
More recent works relating number theory and statistical mechanics can be 
found in Refs. [107-109]. Unlike these authors, we would like to address here 
different type of issues. In particular, we would like to understand better what 
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makes two dimensional models so special from the point of view of number 
theory. Once the answer to this question is found, there is an opportunity to 
think about higher dimensional exactly solvable models using number theory. 

In Section 4.1. two functions G(s) and G^^(s) were introduced in Eq.s (4.10) 
and (4.11). These are the basic functions to which one can relate many two 
dimensional exactly solvable statistical models at criticality. We would like 
to rewrite these functions in the number-theoretic form using simple example 
presented in detail below. 

Consider a ring of complex numbers A := {a + bi \ a, b, e Z}.This is the ring 
of Gaussian integers. Elements of A correspond to points on the lattice in the 
complex plane. For the complex number x — a + bi the norm ||x||is defined as 
usual by 

\\x\\ = a 2 + b 2 = xx. (4.50) 
We define now the (— function for A: 

Since the ring of Gaussian integers has 4 units (that is there are 4 numbers for 
which ||a|| = 1) this should be taken into account when ( A (s) is presented as 
the Euler product 

C A (s) = ^l[(l-\\q\rr 1 (4-52) 

Q 

to be compared with the Hecke L function, Eq.(4.43). The product in the r.h.s. 
of Eq.(4.52) is over prime ideals. To find these ideals, consider the following 
chain of arguments. Suppose a + bi belong to the prime ideal, then a — bi should 
also belong to the prime ideal because if it were possible to decompose a — bi 
into two factors, the same would be possible for a + bi. Hence, a 2 + b 2 has the 
following decomposition into prime factors 

a 2 + b 2 = (a + bi) (a - bi) . (4.53) 

This suggests that either a 2 + b 2 = p or a 2 + b 2 = p 2 where p is the prime 
number. If p = 3 mod 4, then a 2 + b 2 = p cannot be solved, that is p is prime 
element of A. If, however, p = 1 mod 4, then a 2 + b 2 = p is solvable so that 
P = 9i92 with qi = a + bi and q 2 = a — bi. These statements can be trivially 
checked for p = 2. This discussion provides a nice illustration to general results 
stated in Appendix B of our earlier work, Ref.[4]. Using these results we obtain 
the following result for Ca( s )> Ref.[99]: 

<-<•>-< MfO n (t4^) 2 n G^> <««> 

v 7 p=lmod4 v f / p= 3mod4 v y ' 

Here the last term on the right comes from prime elements of A with norm 
p 2 ,the middle comes from two factors both with norm p. Finally, the factor 
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containing 2 comes from the element 1 + i whose norm is 2. Eq.(4.54) can be 
rewritten now as 

C» = 4C(s)L x (s) (4.55) 

with ((s) being the usual Riemann zeta function and L x (s) being the Dirichlet 
L- function 

L xW = £*M, ,4.50) 

n— 1 

where 




0, if n is even 
1, if n = 1 mod 4 
— 1, if n = -1 mod4 



There is an alternative way of presenting the same thing, e.g. see Eq.(2.17). 
Namely, we can write as well 

n=l 

where the factor h{n) denotes the number of ways of representing n as the 
sum of two squares x 2 + y 2 , x, y E Z. Elementary exercise in number theory [99] 
produces the following result for h(n) 14 : 

h(n) = Aj2x(n). (4.58) 

m\n 

Clearly, if instead of the ring of Gaussian integers we would have a ring of 
integers of the imaginary quadratic field, that is a + \^db, we could repeat 
almost word for word all the preceding arguments. This allows us to rewrite 
G(s) of Eq.(4.10) as 

GO,) = E^r (4-59) 

n=l 

where the function h(n) counts the number of integer solutions to the equation 
m 2 + t 2 1 2 = n. Hence, it can be presented as well in the form of zeta function 
analogous to Ca( s )- As for the function G /J .^(s),this is nothing but the special 
case of the function S a (x,y,s) (for y = a = 0) introduced by Kronecker, e.g. 
read Chapter7(section 4) of the famous book by Weil, Ref.[110]. Already Kro- 
necker knew how to connect this function with the theta functions and, since 
these functions can be made to satisfy the KdV equation [72] whose dual is 
related to the Virasoro algebra (e.g see Remark 3.8) it is clear that all results 
of the conformal field theory can be obtained now. One can do much better, 
however. Looking at Eq.s(4.50) and (4.51) one can easily realize that Eq.(4.51) 
is not restricted to numbers of quadratic imaginary field provided that one de- 
fines the norm for the number field. In the case of quadratic field the norm 

14 Hcrc for m = lmod4, x( m ) = b f° r 171 = 3mod4, x( m ) = — b ct c, just like in Eq.(4.56). 
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||x|| = xx. To generalize this result to more general number field it is sufficient 
to define the norm ||x||as the product ||x|| = x\ ■ x 2 ■ ■ ■ x n where X2 r •• i x n are 
"conjugates" of xi.To understand what all this means it is sufficient to realize 
that the algebraic equation 



normally will have n roots. These roots are conjugates of each other. Let Cj 
belong to the field K while the roots to the field extension of K, say E. Let 
<j:E — > E be an automorphism of E which leaves K (that is a ) fixed. If Xiis the 
root of Eq.(4.60), then <t(xi) is also the root of the same equation. The group 
cof automorphisms of E is known as Galois group, Gal(E/K).The Kronecker- 
Weber theorem states that every finite abelian extension of Q (that is finite 
extension of E/Q with Gal(£?/Q) abelian) can be embedded into a cyclotomic 
extension Q(w) with oj being some root of unity. This result underscores the 
importance of the cyclotomic fields discussed in the Appendix A. 

Going back to Eq.(4.51) and replacing A with E we obtain the multidimen- 
sional generalization of G(s). The actual computations will depend crucially on 
the solvability of the Diophantine equation of the type 



Hence, the multidimensional extension of the results of conformal field theories 
(CFT) depends on our ability to solve the Diophantine equations of the type 
shown above. The theory of Diophantine equations is sufficiently developed, 
e.g. see Ref.[99] and references therein. This fact provides hope that multidi- 
mensional extension of CFT can be developed gradually. 

Before leaving this subsection, several comments arc still in place. First, we 
would like to rewrite the Dirichlet L- function, Eq.(4.56), in a different form. To 
this purpose, we need to know that the Dirichlet character x is actually some 
root of unity. If we identify Gal(Q((„)/Q) with (Z/raZ)* (where ( n is root of 
unity and * denotes the multiplicative group of residue classes (such group does 
not contain zero)), then the Dirichlet character modn is the Galois character. 
Since all characters are of the form exp(2iriut /n) = X„(i) it could happen that 
X u (ti) — Xuite) if (^i; n ) — (*2)W) = 1 and t\ = £2 mod/ for some / smaller 
than n.Such / is called conductor of x 15 . It is important to realize that / divides 
n so that x 1S periodic with period /. With this information in our hands, the 
Dirichlet L-function, Eq.(4.56), can be rewritten as [111] 



c + cix + • • • + c n x n = 



(4.60) 



Ml 



E/K 



= n. 



(4.61) 




(4.62) 



with the Hurwitz zeta function being defined as 




< b< 1. 



(4.63) 



Not to be confused with that defined earlier in Eq.(4.29) ! 
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Using this result, Eq.(4.11) can be rewritten with help of Eq.(22) (Chapter 7, 
section 9 of Ref. [110]) in a similar form (which needs some minor modifications 
for higher dimensional lattices) 

G>,i/(s) = H(x, 2s) + H(l - x, 2s) (4.64) 

where x = [i + vr. Using results of Weil's book, Ref. [110], it is not difficult to 
figure out the higher dimensional analogue of Eq.(4.64). In fact, the reader may 
want to consult Refs. [112,113] where this topic is discussed. Their work is to 
be considered also in the next subsection. 

Actually, physically interesting results are only those obtainable as limits 
G'(0) and G'^ „(0) with prime being differentiation with respect to s [82]. Cal- 
culation of these limits leads us directly to the Chowla-Selberg formula for 
periods to be considered in the next subsection. 

4.5 The Chowla-Selberg formula and the Veneziano am- 
plitude 

The Chowla-Selberg formula, Ref. [31], is essentially the same thing as the 1st 
Kronecker's limit formula as discussed, for example, in the books by Weil, 
Ref.[110], and Lang, Ref. [114]. It deals with the unusual reinterpretation of this 
formula due to results of Lerch (Eq.(23) of Chapter 8, section 9 of Ref.[110]). 
Such reinterpretation does not come across easily upon reading Ref. [31]. More 
directly it was obtained in the paper by Ramachandra [75] and was subsequently 
reproduced in Weil's book, Ref. [110]. 

Let K= Q(V— d) be the quadratic imaginary field, i.e. imaginary quadratic 
extension of Q of discriminant -d. Let tk be its ring of integers and suppose 
that tk = Z + Zz with z € K. For such field the Dirichlet function can be 
written as 

L x{s) = j s Y d X{n)H{ T ^,s) (4.65) 

71=1 

to be compared with Eq.(4.62). By analogy with Eq.(4.55) we consider the 
Dedekind zeta function for the field K. Differentiating this function with 
respect to s we obtain: 

R(s) = ^ [C(s)L x (s)} = C'(s)L x (s) + C(s)L' x (s). (4.66) 

Now we have to take the limit s — ► in this equation. Calculation of C'(0) 
is easy and produces ln27r. Calculation of L' x (0) is somewhat more involved. 
Indeed, using Eq.(4.65) we get for L' the following result: 

L' x (s) = [- lnd] L x (a) + ^ £ X (n)H'(^s). (4.67) 

n=l 
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The Lerch formula obtained by Lerch in 1894 ( Eq.(23) on page 60 of Weil's 
book, Ref.[110]) comes to the rescue thus producing 

H'(x,a = 0) =ln^fL (4.68) 

V27T 



The value of L x (0) was known already to Dirichlet: 



M0) = ^31, (4.69) 

where the class number h(\/—d) was introduced in Subsection 4.2 and wq is the 
number of roots of unity in K. Finally, the value of £(0)=— 1/2 while that for 
£'(0) = — 21n27r. Collecting all these results in Eq.(4.66) we obtain: 



R(s = 0) = ~ In 2. + \\ Indl - \ £ X (n) In ^B. (4.70) 

n— 1 



Specializing to those imaginary quadratic fields for which h = 1 the 1st limit 
formula of Kronecker (Eq.(17) on page 75 of Weil's book) is given by 

( K (s^0) = l--A-ln(^(L)) (4.71) 

with.F(L) denoting Siegel-Ramachandra invariant: ^(i) = (27r)~ 24 [iV(L)] 12 |A(r)| 2 , 
Ref.[114],ch-rsl9, 21. Here A(r) = (2tt) 12 g fl (1 — g 2 ™) 24 = (2tt) 12 t?(t) 24 with 

n=l 

t]{t) being the usual Dedekind eta-function while N(L) being the norm of the 
principal ideal. This norm can be easily calculated for the quadratic fields [115]. 
Indeed, if the ideal L=[u>i, a^jthen the different V can be calculated as 



V 



while the same different for the basic lattice L = [1, (r — r)] is given by V —(t — 
r). The norm A^(L) is determined by the rule : N(L) = V/V so that N(L) = w\. 
Keeping in mind that the Dedekind eta function is the modular form of weight 
12 and using Eq.(C.4)one obtains easily: T{L) = (Imr) 12 |A(t)| 2 . 

Differentiating Eq.(4.71) with respect to s and equating this result with 
Eq.(4.70) (for h = 1), we obtain (omiting unimportant constants): 

d 

(Imr) 12 |A| 2 ~ J| T{^) 6wox( - n \ (4.72) 
n=i d 

We would like to compare this result with that known from the CFT. In par- 
ticular, the free boson field theory on the torus is discussed in Ref.[82] , ch-r 
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10, section 10.2. If we let w in Eq.(4.71) to be one then, using this formula we 
obtain at once 

§- s C K (s) |a=o= -2 In [im r ( V (r)) 2 ] = G'(0) (4.73) 

with G'(0) being defined earlier by Eq.(4.10). This result coincides exactly with 
Eq. (10.29) obtained with help of the the path integral for the boson field theory, 
Rcf.[82]. The path integral for such theory is effectively the Euclidean version 
of the bosonic string in the conformal gauge [116]. More details on this topic 
are given in the next subsection. 

In the meantime, use of Eq.(4.72) enables us to provide different interpreta- 
tion of the results known in physics. First, using Eq.s(4.1)-(4.5) and following 
Chowla and Selberg, Ref.[31], we write A(t) = wK 12 with K being one of 
periods of the elliptic curve and w being some (known in principle) algebraic 
number. Taking into account that Imr is also an algebraic number (e.g. read 
Section 4.2) we notice that the combination (Imr) 12 |A| 2 ~ (K) 24 . Therefore, 
using Eq.(4.72) we obtain 

d 

k = cy^n r 0^ x(n) ( 4 - 74 ) 

n=l d 

where C is some constant (algebraic number). This is exactly the Chowla- 
Selberg formula (for h = 1), Eq.(4), page 110 of Rcf.[31]. The factor of 
comes in view of Eq.s(4.1) and (4.70). 

Remark. 4. 8. At this point we would like to remind to our readers that 
exactly the same Dcdckind zeta function (k( s ) used in obtaining the Chowla- 
Sclberg formula was obtained earlier in our work, Rcf.[4], and was interpreted 
as dynamical partition function of 2+1 gravity. 

We would like to demonstrate now how the Veneziano-like amplitudes can 
be obtained with help of such partition function. 

To do so, several things should be kept in mind. First, the number wo is 
equal to 4 for d = 1 (or d=4), it is equal 6 for d = 3 and it is equal to 2 in the 
rest of situations as it was explained earlier, e.g. see discussions after Eq.s(4.18) 
and (4.24). Next, the character x{ n ) m Eq.(4.74) is equal to ±1 since it is the 



Kronecker symbol that is x( n ) = y~ J according to p. 109 of Ref.[31]. With 

this information, let us first look at the simplest case d — 4 having in mind to 
reproduce Eq.(3.37). In this case we need to evaluate the following Kronccker's 
symbols: (|), (§) and (|) (since (|) = by definition). Using Ref.[115], pages 
141-143, we find: (|) = 1, (|) = 0,(|) = —1. In arriving at these results we took 
into account that the discriminant of the field d is not a square- free integer while 
d/A is. Taking into account that ^/tt = T(i) we obtain back the result given in 
Eq.(3.37). Following Chowla and Selberg, consider now another example. This 
time, let d = 7 so that u> = 2. Then, similar calculations produce 



K = Cy/^ 



4 \ "I 2 



r(f)r(f)r(f) 



(4.75) 
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Using Eq.(1.9) it is possible to transform Eq.(4.75) into more comprehensible 
form: 

^=(c-/7r)r(i)r(^)r^). (4.76) 

Looking at the Veneziano amplitude, Eq.(1.12), it is evident that the bracket 
containing sinus factors will not contribute to physically interesting poles. Hence, 
we may rewrite this amplitude as A(a, b, c) ~ r(a)r(6)r(c) keeping in mind that 
a + b + c = 1 (or — 1). In Ref.[117] Gross had shown that it is always possible 
to reexpress K as the product of gamma functions for other d's with h = 1. 
In this sense, the Chowla-Selberg formula, Eq.(4-74), could be considered as the 
Veneziano amplitude. In view of the Remark 4.8., it follows that the parti- 
tion function of 2+1 gravity is capable of producing the Veneziano amplitudes 
and, since this function was not obtained in 26 dimensions, accordingly, our 
Vencziano-likc amplitudc(s) "live" in normal space-time dimensions. This con- 
clusion is in accord with earlier obtained based on Eq.(3.6). 

Obtained result, Eq.(4.74), is inconvenient, however, for use in particle 
physics since it is restricted to a very specific values of parameters, strongly 
depends upon the toroidal geometry and cannot be easily generalized to the Rie- 
mann surfaces of higher genus in spite of several attempts to do so [112, 113]. We 
provided here all these derivations nevertheless because we believe that the 
number-theoretic methods discussed in this work could be potentially useful for 
development of CFT models cither in dimensions higher than two or for models 
defined on Riemann surfaces of higher genus. This is especially important in 
view of the fact that such an extension for higher genus surfaces in mathematics 
literature is only possible for those Riemann surfaces ( and, hence, the Abelian 
varieties, Jacobians, etc.) which admit CM. The theory of CM was developed 
by Taniyama and Shimura [18] and can be found also in the monograph by 
Lang, Rcf. [118]. To our knowledge, the work of Shimura and Taniyama on CM 
had not been used in physics literature so far. 

Before closing this section, we would like to discuss yet another reason for 
discussing the Chowla-Selberg results. 

4.6 Heights of elliptic curves, the Chowla-Selberg formula 
and string theory 

In this subsection we would like to discuss still another connection between the 

Veneziano amplitudes and the number theory. It is based on realization that 

the l.h.s. of Chowla-Selberg formula, Eq.(4.72), can be written in a different 

way. This is possible because the function A(r) is the modular form and, 

hence, following ideas of Eichlcr-Shimura theory discussed in Section 4.3., we 

i 

can construct the invariant differential form a = (A(r)) 12 dz and, with help of 
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such form, the intersection pairing [119]: 



l - j a A a = |A(r)|« dz A dz (4.77) 



C/L C/L 

= |A(t)|® J dxhdy 
A 

= |A(r)|Mmr. 

Following Silverman [89,119], denote the height h(£/Q) of the elliptic curve E 

as 

h(£/Q) = ~hi l - J aha. (4.78) 

C/L 

The above definition is not the most complete. Pedagogically clear explanation 
of the concept of height and its usefulness in arithmetic algebraic geometry is 
given in the review paper by Mazur [120]. Rigorous connection between the 
height and the Chowla-Selberg formula is presented in relatively recent Annals 
of Mathematics paper by Colmez [121]. We deliberately avoid all intricacies of 
this concept 16 since our goal is more modest: we want to connect Eq.(4.78) with 
string theory. 

Following Ref.[8], the Weil-Petersson fundamental (1,1) Kahler form ojw-p 
on the Tcichmiiller space of the torus (and/or the punctured torus [9]) is given 
by 

u W -p = j dl ~ A df ■ (4.79) 

4 (Imr) 

At the same time, the genus one partition function for the bosonic string in 26 
space-time dimensions is given by [123,124]: 

Zi = const [ UW ~ P (4.80) 

! (Imr) 12 |A| 2 



M 



i.i 



const j u>w-p exp(24/i(£/Q)) 



Mi 

with integration domain taken over the moduli space A^i^of once punctured 
torus. In our previous work, Ref.[64], we had discussed extensively the above 
integral for the case when the height function is equal to zero. This leads imme- 
diately to the Kontsevich-Witten matrix models, etc. Since the height function 
is closely connected with the Arakelov theory [119,121], naturally, extension 
of Eq.(4.80) to the Riemann surfaces of higher genus involves elements of this 
theory. This indeed was accomplished by Manin [125] and Bost and Jolicoeur 
[126] ( see also paper by Smit [127] ). The above formula and its higher genus 

16 These can be found easily in the literature just cited, e.g. sec Rcf.[122] in addition. 
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generalizations contains several important drawbacks. First, the modular form 
A(t) is a cusp form. That is, it can exist only if the torus is cusped, i.e. it has 
at least one puncture. (The dimension (and, hence, the number) of cusp forms 
surely depends upon the actual number of cusps). This topology is the sim- 
plest possible for the case of 2+1 gravity as it was extensively discussed in our 
earlier work, Ref.[4]. But, in string theory the puncture is normally associated 
with some particle via vertex operator insertion. By definition, the string 
partition function without external sources (or sinks) should not contain 
punctures. Second, the expression for the height makes sense in the Diophan- 
tine geometry dealing with finite number of points on the algebraic curves. The 
Chowla-Selberg formula, Eq.(4.72), is in perfect agreement with such statement 
since it involves discrete t's coming from the imaginary quadratic fields. Use 
of Eq. (4.80) destroys the complex multiplication option and, hence, disconnects 
the l.h.s. of Eq.(4.72) from the r.h.s. 17 If, however, the complex multiplication 
is destroyed (not used, ignored) then, according to Weil, page 40, paragraph 7, 
Ref.[110 ], the very basic transformation law, Eq.(C.7), of Appendix C for the 
modular functions is violated. Third, the genus zero partition function seem- 
ingly does not require integration over the moduli space [7]. This, however, 
is suspicious. Indeed, the open string world sheet for genus zero case is repre- 
sented by the upper Poincare halfplane H which is a model of hyperbolic (that is 
non Euclidean) space. The minimal number of open external states is four and 
they appear as insertions at the boundary of H (e.g. sec Fig. 1.15(c) of Rcf.[7]). 
Technically speaking, these are just the cusp points. The four times punctured 
sphere is a hyperbolic surface whose cover is H with four cusps. The locations 
of cusps is not fixed. The motion of the cusp points is described by the motion 
in the parameter, i.e. in the Teichmuller, space. Nag had demonstrated [9] that 
the Teichmuller space of the four times punctured sphere is the same as of 
once punctured torus. Notice, that in view of these arguments there is no need 
to discuss separately open and closed strings. This explains why the Veneziano 
amplitude makes more sense than the Shapiro- Virasoro. If one accepts the 
above arguments, then, 2+1 gravity and string theory are in agreement with 
each other so that formally Eq.(4.80) is a partition function for "time-reduced" 
2+1 gravity. This was noticed already in our earlier work, Ref.[64] and further 
developed by Krasnov [131] 18 . Fortunately, there is much more efficient alterna- 
tive route of obtaining the multiparticle Veneziano amplitudes. It is presented 
in the next section. 

17 This can be repaired if the traditional Tcichmullcr/moduli space is replaced by its p-adic 
analogue [128]. Even then, as the recent developments indicate [129,130], one still will end up 
with the p-adic version of the results presented in Section 5 below. 

18 For more recent references, please, see the footnote # 3. 
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5 Fermat's hypersurfaces and Veneziano ampli- 
tudes 



5.1 The Picard-Fuchs equations 

In Section 3 we had shown that the standard Veneziano four particle amplitude is 
just one of the periods of the Jacobian variety for the Fermat curve. For physical 
applications it is of interest to obtain the multiparticle amplitudes. This can be 
done in two ways: cither by explicit calculation based on generalization of the 
results presented in Section 3 or by considering some sort of equations whose 
solutions produce the desired periods. The situation in the present case is very 
similar to that encountered in mirror symmetry calculations [57]. Because of 
this similarity we will be brief in discussing the second option. But, we feel, 
that this option should be mentioned explicitly in the text because it provides 
some additional insight into physical and mathematical meaning of periods and 
their calculation. 

We begin with the simplest example borrowed from the pedagogically written 
paper by Griffiths [132]. He considers calculation of the period of the following 
integral along the closed contour T in the complex z-plane: 

r 

Since this integral depends upon parameter A the period 7r(A) is some function 

of A which can be determined as follows. Differentiate 7r(A) in Eq.(5.1) with 

respect to A. This produces: 7r'(A) = § jr0^p ■ The combination 

r 

A7r'(A) + 7r(A) =0 (5.2) 

produces the desired differential equation which enables us to calculate 7r(A).This 
simple result can be vastly generalized to cover the case of periods of integrals 
of the type 

= / P J" / Zl -'" ,Zn \ dz 1 A dz 2 ... A dz n (5.3) 
/ Q(z lt ...,z n ) 
r 

The equation Q(z\ t ...,z n ) = determines the variety. It may contain parameter 
(or parameters) A so that the polar locus of values of z's satisfying equation 
Q = depends upon this parameter. By analogy with Eq.(5.2) one obtains 

k d n 

^P„(A)- 5r O(A) = 0. (5.4) 

n=0 

This is the Picard-Fuchs equation for periods. To solve this equation one needs 
to know the explicit form of polynomials P„(A) . In general this is not an easy 
task as it was demonstrated by Manin [133] many years ago. And, because this 
is not an easy task, we believe, that calculation of periods using generalization 
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of section 3 is more efficient. Nevertheless, method of differential equations 
teaches us things which are not immediately apparent when direct calculations 
of periods is made. 

In Section 4.1. we noticed that for the elliptic curves the periods can be 
calculated as solutions of the hypergeometric equation, Eq.(4.6). Surely, this 
equation is of Picard-Fuchs type. However, in Eq.(4.6) there is no explicit 
differentiation with respect to parameter A and, hence, now we would like to 
correct this deficiency. Fortunately, this task was accomplished by Manin [133]. 
We follow Ref.[134], however, which is more elementary. The Legendre form, 
Ref.[83], page 179, of the elliptic curve £ is given by y 2 — x(x — l)(x — A). The 
invariant differential form lo on £ is given by 

/ \ \ dx dx 
Lo(X) = — = — . (5.5) 
2y 2y/x(x-l)(x-\) 

Differentiation with respect to A produces 



A(A - 1) J^(A) + (2A - 1) A W(A ) + = d{ % (5-6) 

i i 

Integrating this relation along the contours 2 J = <f a and/or 2 J = <L we obtain 

A a 

the equation for periods 



A(A- l)^n(A) + (2A-l)An(A) + ln(A) = (5.7) 



in accord with Manin, Ref.[133]. Let us recall [83] that for the elliptic curve in 
the Legendre form the j— invariant is known to be 



, 8 ( A 2 -A + 1) 3 
A 2 (A-1) 2 



Use of this invariant allows us to compare elliptic curves between each other. 
In particular, if two curves are isogenous they have the same j— invariant. This 
means that j— invariant (or parameter A) describes the family of all elliptic 
curves and, hence, A can be identified with the point in the moduli space of el- 
liptic curves. This means that the differential operator in Eq.(5.7) is defined on 
the moduli space of elliptic curves and it effectively describes the deformation of 
complex structure for such curves. Since the moduli space is non Euclidean, the 
differential operator in Eq.(5.7) is actually defined in some "curved" space char- 
acterized by some connection (Gauss-Manin connection). The commutator of 
covariant derivatives on such space produces zero curvature (so that connection 
is flat but not trivial nevertheless). Surely all this machinery can be extended 
to more sophisticated algebraic varieties. A good summary of results applicable 
to general case can be found in Refs. [36,134,135]. In this paper we are not go- 
ing to use these results however and, hence, the above discussion is only meant 
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to indicate the alternative route to reach the same destination. The necessity 
to mention about such a route comes from the fact that it actually provides 
the unified description of both conformal and "string" theories as men- 
tioned briefly already in Section 4.1. Although the Picard-Fuchs equation for the 
B-P type singularities can be easily derived (as discussed in the Appendix D), 
this route for obtaining the multiparticle Veneziano amplitudes is not unique. In 
the next subsection we discuss another route closer in spirit to that envisioned 
by Landau [24] and further developed by Pham [25] and others. 

Conjecture: Based on the results of Pham [25], Lefshetz [136] and Griffiths 
[132], and those presented in this work, it is natural to expect that any scat- 
tering amplitude (correlation function) in particle physics (in conformal field 
theories) should come out as solution of some kind of the Picard-Fuchs equa- 
tion and, hence, should be interpreted as one of the periods of differential forms 
" living" on the moduli space appropriate for the particular scattering problem. 
In the case of conformal field theories this conjecture could be considered as 
actually proven in view of the work of Schehtman and Varchenko [137] (to be 
discussed briefly below in connection with hyprplane arrangements) while in the 
case of particle physics this work along with those on mirror symmetry [57] (and 
references therein) could be considered as first steps towards its proof. 

5.2 Veneziano amplitudes from Fermat hyper surfaces 

To obtain the multiparticle Veneziano amplitude we have to generalize results 
of Section 3.2. It is appropriate to mention at this point that development of 
the multiparticle Veneziano amplitudes leading to string theories (as discussed 
in detail in Ref. [19]) differs markedly from presentation which follows. Our 
presentation is based on works of Deligne [37] and Milnor [26]. Additional 
mathematical details and references can be found in Ref. [38]. 

In section 3.2 the main object of interest was the differential 1-form, Eq.(3.24), 
leading to the Veneziano amplitude, Eq.(3.33). Now we need to find the analo- 
gous form for the Fermat hypersurface, Eq.(3.15), with ao = • • • = a n = N. To 
do so, we would like to remind our readers about some facts from the theory of 
Riemann surfaces to be generalized to more complicated case of Fermat hyper- 
surfaces. In the case of Riemann surface R there are three types of differential 
forms lo. For the first and second type the period integrals J to are non-zero 

7 

only for non-null homologous cycles 7 on J? whereas for the third type one 
obtains the so called residual periods associated with poles of differential forms 
oj.They can occur for null homologous cycles as well and, hence, are of no in- 
terest to us [138]. Suppose now that we have a p-form to on a manifold X. 
If we integrate such a form over p-cycle, then one obtains a period (usually a 
complex number) depending only on a homology class of the cycle provided that 
lu is closed (that is dui — 0).The homology classes form p th homology group of 
X with complex coefficients, H p (p,C).The dual to it is p th cohomology group 
H p (p,C)~H p (p,C)*. A closed p-form to yields a cohomology class [to] €H P (X,C). 
The Stoke's theorem shows that this class remains unchanged if we add to it 
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an exact form, that is differential dtp of the p — 1 form ^.This circumstance al- 
lows us to introduce the equivalence relation and, hence, to study the quotients, 
e.g. the de Rahm p th cohomology group H^(X) = C P (X)/EP(X), with C p 
and E p being respectively the vector spaces of closed and exact p-forms with 
complex coefficients. According to de Rham, Ref.[138,139], there is a canoni- 
cal homomorphism: H^(X) — >H P (X,C). Consider first how all this applies to 
the Riemann surface R which is complex one-dimensional manifold X. For such 
manifold H^^X) is isomorphic to the vector space of holomorphic 1-forms. 
Each holomorphic 1-form ui is closed. This is so because dui is a holomorphic 
2-form which is equal to zero automatically since X is one dimensional. From 
topology it is known that for the Riemann surface of genus g the cohomology 
group H 1 (X,C) has dimension 2g . Meanwhile the dimension of holomorphic 
1-forms is g so that the holomorphic 1-forms yield only 1/2 of cohomology. 
To repair this situation one needs to consider not only the cohomology classes 
coming from the holomorphic forms but also coming from the anti holomorphic 
forms. For a general manifold X, following Hodge, one introduces the differential 
(p, q)- form lo (in local complex coordinates z\,...,z n ) such that 

u> =^2a h ...^^ ... ijq dzi A • • • A dz ip A dz h A • • -Adz jq (5.9) 

Each m-form on X can be expressed uniquely as a sum of (p, g)-forms with 
p + q = m. 

The following theorem plays a very important role. 

Theorem 5.1. If X is compact complex manifold, then it is Hodge manifold 
( or Hodge-type manifold) if and only if it is a projective algebraic manifold at 
the same time. 

This statement is known as the Kodaira's embedding theorem, Ref.[139], 
ch-r 6. Evidently, the projective algebraic manifold must be complex. 

Definition 5.2. A complex projective manifold X is a nonsingular sub- 
variety in CP" consisting of the common zeroes of a system of homogenous 
polynomial equations in coordinates (z , z n ) of C n+1 . 

On C n+1 one can introduce the Hermitean scalar product (which is just the 
finite dimensional analog of what is being used for the Hilbert space of quantum 
mechanics) : 

n 

< z,z >~^^ZiZi. (5.10) 

i=0 

With help of such defined scalar product it is possible to recover the metric and 
the curvature for such spaces. Moreover, with little efforts this can be extended 
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to CP™ [139]. If the Hermitian metric h of the projective manifold X embedded 
in CP™ is given by 

h = ^ h llv (z)dz li ® (5.11) 
then, the fundamental (1,1) form Q is given by 

9.= ^Y^h^{z)dz^Adz^. (5.12) 

Definition 5.3. The Hermitian metric ft, on X is of Kahler type if £1 is 
closed, i.e. dQ = 0. 

Not all complex manifolds admit the Kahler metric 19 . Fortunately, those 
which admit have direct physical relevance. In particular, all known physi- 
cally meaningful symplectic manifolds of exactly integrable finite dimensional 
classical dynamical system admit Kahler type metric and Q, form. [40,141,142]. 
Using Kodaira's embedding theorem we conclude that all physically interest- 
ing projective manifolds of Kahler type are also manifolds of Hodge 
type 20 . 

In particular, let us consider the case of Riemann surface of genus g. It can 
be mapped into g dimensional complex torus (the Jacobian) under appropriate 
conditions (e.g. it must be a polarized Abelian variety) [118,134,139,140]. And 
such torus surely admits Kahler metric and, hence it is a Hodge-type manifold 
with Hodge decomposition: 

H 1 (X,C)=H 1 -°(X)e ^(X). (5.13) 



Remark 5.4. The condition that an abelian variety is of complex multipli- 
cation type can be interpreted as a choice of basis for the above decomposition, 
Eq.(5.13),[143]. 

Remark 5. 5. The situation becomes more complicated if complex hypersur- 
face possess some kind of singularities, e.g B-P type discussed in Section 3.2. 
In this case the mixed Hodge structures should be considered instead of pure 
Hodge structures as shown by Deligne [37]. Systematic up to date exposition 
of this topic can be found in Ref.[135]. Mathematically rigorous and physically 
accessible treatment of these structures can be found in excellent review paper 
by Varchenko [144] (Appendix D). 

Complex multiplication (discussed in Section 4.2) is closely associated with 
manifolds of Hodge type. Fortunately, this is indeed the case for both the 
Format curves considered in Section 3 [34] and for the Fermat hypersurfaces [37] 
to be considered below in this section. Construction of Hodge type manifolds 
admitting complex multiplication is concisely described in the paper by Weil 
[33] and also in his book, Ref.[140]. More recent and much more advanced 

19 Comprehensive and readable account of Kahler manifolds is given in Andre Weil's book, 
Ref.[140 ]. The book by Wells, Rcf.[139], is also an excellent source. 
20 For additional details, please read Section 6. 
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presentation of this topic involving theories of motives, Tannakian categories, 
crystalline cohomology, etc. can be found in Refs. [38,145]. 



Consider now the projective form of Eq.(3.15) adopted to the Fermat case, 

z? + ■■■ + *£ -^=0. (5.14) 

Let us consider it as a variety over the cyclotomic number field K=Q(cxp{2iri/N}). 
Then, using results of Appendix A we can rewrite it as a hyperplane in PC™, i.e. 

zq H Vz n - z n+x = 0. (5.15) 

Such a hyperplane can be presented as a point in the Grassmanian as we had 
discussed in our earlier work, Ref.[64]. Each hyperplane, Eq.(5.15), is just a 
result of a particular solution of the associated system of cyclotomic equations 
(e.g. see Eq.(A.3)). Naturally, there are many other solutions and each of them 
will be represented by its own point in the Grassmanian. One can think about 
the transformations which connect different hypcrplancs. Such transformations 
are presented by the matrices G belonging to the rotation group. Among these 
transformations might be those, call them N, which have the fixed point. If this 
is the case, then, one should consider a qotient G/N. This quotient is essen- 
tially what is known in the literature as the Griffiths domain [139, 146, 147]. The 
Griffiths domain is the moduli space for periods. We shall demonstrate 
momentarily that each hyperplane is associated with its period (that is with 
the multiparticle Veneziano amplitude) so that the totality of hyperplanes are 
in one to one correspondence with the totality of periods. This fact is discussed 
further in Section 5.3.2. and Appendix D. In principle, if one can construct the 
matrix of periods then, one can determine the Griffiths domain as well. To keep 
focus of attention of our readers on physics, mathematically relevant issues of 
physical importance are discussed in the Appendix D. This allows us to come to 
the main purpose of this work- calculation of multiparticle Veneziano amlitudes. 

To this purpose, we need to use the Corollary 2.11. from Griffiths work, 
Ref.[132]. It is related to the integral, Eq.(5.3).To explain things better, we 
adopt the explanation of this Corollary using discussion from the Brieskorn's 
book, Ref.[138], page 647. 

Let (xo, x n ) be homogenous coordinates of the point in the projective 
space and ( z\,...,z n ) be the associated with them coordinates in the affine 
space with Zi — Xi/xo , then the rational n-form is given by 

gjfWgnK t Adz2 AdZn / 5>16 n 

q{z lt ...,z n ) 

with rational function p(z)/q(z) being a quotient of two homogenous polynomi- 
als of the same degree. Under substitution — Xj/x , the form dziAdz 2 ---/\dz n 
changes to 

n 

dzi A dz<i... A dz n — x < - n+1) Y^(— l) 4 Xidx A ... A dxi A ...dx n (5-17) 
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where the hat on top of Xi means that it is excluded from the product. Define 
now the form 

n 

ujq := ^(-l) 4 Xidx A ... A dobi A ...dx n . (5.18) 

i=0 

In order to rewrite the form, Eq.(5.16), in terms of x coordinates it is sufficient 
to multiply the denominator of p(z)/q(z) by x^ +l and to replace z by x in the 
quotient. Thus one obtains a homogenous expression 

in which the (degree of Q) = (degrec of P) +(n+l). This is Griffiths result, 
Corollary 2. 11. of Ref.[132]. Conversely, each homogenous form, Eq.(5.19), yields 
an affine rational form, Eq.(5.16), by substitution: x = l,Xi = Zi. 

Next, consider the set X(S 1 ) of numbers which we call the spectrum of the 
critical point (Appendix D) of the B-P singularity, Eq. (3.7). It is given by 

XiS 1 ) = {ae (Z/NZ) n+2 = (Z/NZ) x • • • x (Z/NZ) | a=(a , a n+1 ), £). en = OmodTV} 

(5.20a) 

Let < cii > denote the representative of a,i in Z such that 1 << >< N— l(e.g. 
see Eq.(3.17a)). Also, define the average 

<a>=-^E 2 <a*> (5.20b) 

then, for the Fermat hypersurface the differential form ui, Eq.(5.19), can be 
written as 

„<ao>-l <a n +l>-l 
« + X n+l) 

with t^o defined by Eq.(5.18) (with n replaced by n + 1). With such form at our 
disposal it will be advantageous for us to use the affine representation of this form 
using results just described. To this purpose, by analogy with developments in 
Section 3.2, e.g. see the discussion after Eq.(3.15), we make change of variables 

: Zi = tf cxp(^p) = tfC , provided that = 1,U > 0, thus forming the 

n + 1 simplex A. Such change of variables allows us to consider a map /(A): 
A — > V(F) from the simplex A to the affine variety V(F) 

Y N + ■■■+¥» = 1, Y t = Xi /x n+1 = m. (5.21) 

This map can be made in many ways: for example we can suppress all £'s , then, 
suppress all, except one, then, all, except two, etc. and still will get the desired 
mapping. To avoid gucssings, we need to reobtain back the result, Eq.(3.31), for 
the simplest four particle Veneziano amplitude. Clearly, such amplitude should 
be obtainable from general result, Eq.(5.20), i.e. it should come up as period of 
the integral J co. Such an integral contains a pole. There should be a procedure 
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generalizing that known in the standard complex analysis of one variable of 
calculating residues of multidimensional complex integrals. Fortunately such 
procedure exist and was developed by Leray, Pham and others [136,148] in 
connection with Landau's work on analytical properties of scattering amplitudes 
of Feynman's diagrams [24]. With such information at out disposal, we would like 
to reobtain the simplest Veneziano amplitude now. Although results of Deligne 
[37] are of great help in accomplishing this task, unfortunately, we cannot borrow 
them entirely for (physical) reasons which will become clear upon comparing of 
our derivations (below) with those by Deligne. 

We begin with Eq.(3.31) written (up to constant factor of N) as follows: 

i = ^C^(i - C')(i - C s )(i - C r )r(a)r(&)r(i - a- b) (5.22) 

-t 

This form clearly calls for elimination of the factor £ 2 for reasons of symmetry. 
This was actually done in Eq.(3.33). This time, we would like to do such elimi- 
nation more systematically. Hence, for the time being we shall keep this factor 
dropping it at the end. Eq.(5.22) can be presented equivalently in the following 
form: 



T(a + b)I = (1 - C)(l - C) J J dx 1 dx 2 x a 1 - 1 x b 2 - 1 exp{-x 1 -x 2 ) (5.23) 



00 1 
= (1 - 0(1 - C) J dtt a+b - 1 exp(-t) J diiaj-^l - x^- 1 


In going from the first line to the second we have introduced new variables : 
X\ = X\t , x 2 = x?t , x\ + X2 = t implying that X\ + £2 = 1 and the Jacobian 
of transformation equal to one. The above can be rewritten therefore as 

1 

/ = (1 - C s )(l - C) J dxtotrHl - ^f- 1 (5.24) 


Looking at Eq.s (5.17) and (5.20) and comparing with Eq.(5.24) we can formally 
write: 7 = (1 — C s )(l — O J 7 w • And, indeed, this expression coincides with 
Eq.(3.31) in view of Eq. (3. 28). These results can be vastly generalized now. First, 
we would like to notice that if we would make a change: r — * — r , s — > — s, the 
arguments leading to Eq.(5.24) will remain unchanged. Therefore, this means 
that this result can be used for calculation of physical Veneziano amplitude, 
e.g. see Eq.s. (1.6) and (3.33). Second, based on this observation, we can also 
change £ s into (~ s and £ r into without changing a and b in the integral 
of Eq.(5.24). Vice versa : we can change a and b without change of phase 
factors. To illustrate why this is possible consider as before change of variables: 

Zi = t™ exp(^). In the integral, Eq.(5.24), we choose z\ =x\ = ((where ( 
means that the £ factor may or may not be present. It may be present if we 
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bring the external phase factors inside the integral, Eq.(5.24)) and, accordingly, 
j_ . 

Z2 = X2 = t^C so that for the Fermat curve in the affine form: + z 2 = 1, 
we obtain back equation for the simplex: t\ + t 2 — 1. Hence, we may consider 
the most general case of the integral of the Veneziano-type: 

Iv { F)= /n ai --'0 1 f i A---A^. (5.25) 

A 

Naturally, it is determined up to a constant. The constant is easily determined 
according to the following rules: 

a) if the number of phase factors is even the overall sign of the integral is 
"+", otherwise, it is 

b) the total phase factor is just £ ai + a2H ^ a n+i ; wnere <jj means that this 
factor may be either zero or some rational fraction; 

c) the sum J27=o ai ~ ® moc ^ ^ > 

d) upon change of variables:! 7 ; = (t~ , the overall c-factor c = (-^-)™ 
emerges (and can be dropped); 

f) the resulting integral 



I v(F)=J V • • • C+ + i — A • • • A — 



with hi — di/N can be calculated in the spirit of ingenious trick by Deligne 
[37] which, of course, can be as well inferred directly from Eq.(5.23).To this 
purpose, we multiply both sides of the last equation by r(^"j" 1 bi)I v ^ — 
r(Er=/ b i) ft h i--- tn+i 17 A • • • A . For the r.h.s. we get 

A 

"+ 1 • - f 6j .6 n+1 dh dt n 



oo 



dt n 

A • • • A — — 

trt. 



A 

Let now Sj = tti so that, as before, t = J27=i s i> s i — O.This then produces the 
following integral 

oo oo 

[ ... /* CXD '_ V" +1 S ■) s bl c b " +1 dSl A A dSn A dSn+1 
J J cx p( 2. i=1 s *> S i Sl s n s n+1 



= T(h) ■ ■ ■ T(b n+1 ) 
so that finally we obtain the contribution to the Veneziano amplitude 

n+1 

n r(6i) 
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The total amplitude is the sum of the above with the appropriate phase factors. 

It should be clear by now why it is permissible to multiply both sides of 
Eq.(3.33) by factor of £ 5 (and, of course, to do the same in the most general 
case). 

Remark 5.6. Eq.(5.27) can be found in Gross paper, Ref.[34], where it is 
given without derivation. The role of complex multiplication in obtaining this 
result is strongly emphasized in his paper. The same equation was formally 
obtained by Edwards [29] already in 1922. 

Remark 5.7 In accord with Section 3.2, different combinations of a^s corre- 
spond to different periods of Fermat hypersurface. Each such period corresponds 
to a point in the Grassmanian so that different periods correspond to different 
points. This makes calculation of the Veneziano amplitudes similar to our 
earlier calculations of the Witten-Kontsevich averages, Ref.[64]. 

Remark 5.8. There is yet another method of obtaining the results of this 
subsection. It is based on the asymptotic analysis of complex oscillatory in- 
tegrals as discussed in the book by Arnol'd, Gusein-Zade and Varchenko, Rcf. 
[36]. Fortunately, the results of this reference not only useful for reobtaining 
the multiparticle Veneziano amplitudes (as discussed in the Appendix D) but 
also they provide a very economical way of getting all the results associated 
with the mixed Hodge structures. These results will be utilized in Section 5.3.2. 
in connection with the Hodge spectrum which is in one-to-one correspondence 
with the particle mass spectrum. 

5.3 Dynamics of Seifert fibered phase and Veneziano am- 
plitudes 

5.3.1 Connections between zeta function and Alexander polynomial 

In our previous work, Ref. [4] , we had considered dynamics of the Seifert fibered 
phase of 2+1 gravity (e.g. see Appendix 3 for definitions) associated with so- 
lution of the Witten-Kontsevich model. In this work we have reobtained the 
partition function (the p-adic version of this function) of 2+1 gravity in Section 
3.1., Eq.(3.6), working with dynamics of gravity in this phase. Earlier, we ar- 
gued that such partition function is associated also with the p-adic analogue of 
the four particle Veneziano amplitude, Eq.s(l. 15), (1.16b), (1.17). In Section 3.2. 
we introduced the generalized Alexander polynomial, Eq.(3.20). According to 
Milnor [26] and Brieskorn [27], such polynomial describes the multidimensional 
fibered knots and links. Since according to the Corollary 3. 7. the spectrum of 
particle masses emerging as poles in four particle scattering amplitude is in one- 
to- one correspondence with the appropriate Alexander polynomial describing 
fibered knots and/or links in 3 dimensions, it is natural to expect that the same 
will be true for the multiparticle scattering 21 . 

21 As the development presented below indicates, this option is less convenient than that 
emerging from the correspondence of mass spectrum with the Hodge spectrum. 
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Remark 5.9. It should be noted at this point that the Alexander polyno- 
mial may describe more than one knot/link. This happens in the case of cable 
knots/links, Ref.[149], pages 121-122. In such cases one might think of one-to- 
one correspondence between the mass spectrum and knot/link cable families. 
This option is not going to be explored in this paper in view of the remarks 
made in the footnote below. 

To proceed, we need to explain why two apparently different equations, 
Eq.(3.13) and (3.18), describe the same Alexander polynomial for the Seifert- 
fibcred knots/links. The connection between the above equations could be 
thought as a corollary to a very deep result of Grothendieck [62] which can 
be formulated in the form of the following 

Theorem 5.10. The Alexander polynomial for algebraic knots/links em- 
bedded into S 2n+1 is equal to the product of cyclotimic polynomials. 

Remark 5. 11. For the purposes of this work it is sufficient to keep in mind 
that the link K=V B -p(f) n S 2n+1 introduced in Section 3.2. is algebraic [149]. 

Remark 5.12. The Alexander polynomial Ajc(t), Eq.(3.18), can be ob- 
tained from the ratio of cyclotomic polynomials given below 

-!)(*-!) 

Ak(<) - (tp-w«-iy (5 ' 27) 

Indeed, using equations (A.6)-(A.8) of Appendix A we obtain for instance, 

^- 1 = (*- ix**- 1 +--o = ![(*-&). 

i 

with £ 4 being one of the roots of unity (including one). Use of such type of 
products in Eq.(5.27) produces, indeed, the desired result, Eq.(3.18). For 3- 
dimcnsional Seifert-fibered knots /links the result, Eq.(5.27), had been obtained 
by Le Trang [150] and, independently, by Sumners and Woods [151]. The re- 
sult, Eq.(5.27), leads to the following additional observation originally made by 
Milnor [26]. It appears, that with help of Eq.s (2. 30), (3. 16), (3.17) one can 
connect the Alexander polynomial and the Weil-type zeta function. Earlier, in 
Eq.(3.6) we had obtained such zeta function for the p-adic version of Veneziano 
amplitude. Moreover, for the torus the Weil- type zeta function, Eq.(3.35), was 
obtained earlier by Dwork [102], e.g. see Remark. 4. 7. Looking at Eq.(3.35) it is 
educational to rewrite it in the following form: 

1 _ («-Dfe«-l) (528) 



Zk{u) (au — l)(qu/a — 1) ' 

Assume now that q—t mn ,a = t' m so that q/a = t n . In addition, assume that 
u = u n { or u = u m ,etc). Then, using cyclotomic-type expansion ( dispayed 
above) and letting u = 1 at the end of calculations produces back the Alexander 
polynomial, Eq.(3.18), for the torus knots. This then implies that, at least for 
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torus knots, ZK(t)Ajt(t) = 1. This result is not totally surprising if one recalls 
[152] that for a well behaving matrix A one has the following identity 

00 j.n 

[det(I - tA)}' 1 = exp{J2 ~ (trA n )} (5.29) 

n=l 

Looking at Eq.(3.13) and recalling that dct M=±l (e.g. read Section 3 of 
Ref.[4]) and comparing the r.h.s. of the above identity with Eq.(2.25) one 
obtains the desired relationship between the zeta function and the Alexander 
polynomial. In Section 3 of Ref . [4] we have demonstrated that the monodromy 
matrix M is indeed involved in surface dynamics. 

Remark 5.13. The cyclotomicity (Theorem 5.9) implies that such dynamics 
is characteristic only for the Seifert fibered phase (regime). This is in accord 
with results of Milnor's book [26], e.g read chapter 10. 

The question remains: will such relationship hold for surfaces more com- 
plicated than torus? Although the answer to this question can be found in 
Milnor's book [26], we prefer to approach it differently. To this purpose let us 
notice that the arguments used in obtaining of Eq.(5.28) can be also applied to 
zeta function, Eq.(3.6), the p-adic Veneziano amplitude. In this case we obtain 

ZJW)= 1-7= n (*-0 = Ako(*), (5.30) 

where Ak (t) is the Alexander polynomial for the " empty" knot [58] .This result 
follows from Milnor's general theory discussed in Section 3.2. Indeed, let in 
Eq.(3.7) there is only one term , i.e./(z) = z".We can still construct a circle 
mapping, Eq.(3.9), and after applying the same type of arguments leading to 
the Alexander polynomial, Eq.(3.20), we will get Ak (t). Just made observation 
can be vastly generalized with help of the fundamental theorem of Thorn and 
Sebastiani [153]. The arguments by presented in Ref. [154] are based on this 
theorem and summarized in the following additional theorem (Oka, Ref. [154]): 

Theorem 5.14. Let f be a polynomial in C n xC m such that f(z,w)=g(z)+h(w) 
for each (z,w)& C™ x C m where g(z) and h(z) are weighted homogenous polyno- 
mials m C" and C m respectively. Let X = / _1 (1) C C" x C m ,y = 5 _1 (1) C 
C" and Z=KT X (1) C C m . Then there is a natural homotopy equivalence a be- 
tween X and Y*Z where Y-kZ is the join of Y and z with strong topology. 

As a corollary, the same author obtains the following result: 

Theorem 5.15. Let g, h and f be the same as in Theorem 5.14- and assume 
that g and h have isolated critical points at the origin. Then the characteris- 
tic (Alexander) polynomials of the associated fiberings A g (t), Ah(t) and Af(t) 
satisfy the equation 

A / (t) = A fl (t)*A h (t). (5.31) 
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That is letting A g (t) = Yl^t — Aj) and Ah(t) = fl^t — /liJ we should obtain 

A/(i)=n ilj (*-A i M J )- 

Remark 5.16. The concept of a join could be found in any book on topology, 

e.g. see Ref. [60]. We actually have effectively used it already when we had 
considered the map A — > V(F) as discussed before Eq.(5.21). Rougly speaking, 
the concept of a join could be understood based on the following example taken 
from Ref.[155] (see also Pham[25] and Milnor[26]). Let Z/p and Z/q be finite 
cyclic groups consisting of all p-th, respectively q-th roots of unity. The join 
J=Z/p*Z/q could be thought as totality of all vectors (s£,tr]) € C 2 with s,t 
> 0, s + 1 = l,and £ €Z/p, r\ £Z/q . For the circle map of the type 

, , , zP + 
ip{z,w) 



zP + wi\ 

it could be shown [26], that J is the deformation retract of the fiber ip~ (1). 

Remark 5. 17. From Theorems 5.14. and 5.15 it follows that the Alexander 
polynomial of the multidimensional fibered knot / link can be constructed with 
help of the Alexander polynomials for empty knots, e.g. A g (t) = Yli(t ~ A^) , 
Eq.(5.30), associated with zeta functions of circle maps related to periodic sur- 
face automorphisms characteristic for the Sefert-fibered phase (Appendix B) 22 . 

With all useful information just discussed we still have not provided an 
answer to the question we had posed earlier about the relationship between 
the Alexander polynomial and zeta function for multicomponent links. Fortu- 
nately, such problem was discussed by A'Campo [156] so that it remains only 
to illustrate his findings. In particular, for the B-P type singularity, e.g. see 
Eqs.(3.7),(3.8), the Alexander polynomial A(i) should be 

A(t) = [(t - l)" 1 ^ - l)x( s «)](- 1 )" (5.32) 

with x(Sn) being the Euler characteristic of the complement of the Fermat hy- 
persurfacc in CP™. A'Campo provides the following result for the characteristic 

: x(Sn) — 1 ~^ 1 ~/y V ^ • Next, if the Milnor number (that is the multiplicity of 

the critical point) is given by 

H = (-1)"[-1 + N X (S N )] = (N- l)"+\ (5.33) 

then the zeta function Z(t) and the Alexander polynomial A(t) are connected 
through 

A(t)=f[^Z(l)](- 1 ) n+1 . (5.34) 

Consider now a special case : N — 3, n — 2, treated by A'Campo. It is relevant 
for 4-particle Veneziano amplitude discussed in Section 3.2. The treatment of 
general case will become obvious upon completion of discussion for this special 

22 Essentially the same conclusions can be reached using the laguage of singularity theory 
(Appendix D), e.g. see the Fubini's theorem, Eq.(D.7), and the discussion which follows. 
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case. The Alexander polynomial for this special case can be written at once as 
follows: 



here p is the same as in Eq.(3.19). The reader familiar with knot theory might 
recognize immediately in the above expression the Hosokawa-type Alexander 
polynomial for 2 component link. To prove that the above expression is, indeed, 
valid for a link, it is sufficient to invoke the theorem by Murasugi, Ref.[157], 
page 117: 

Theorem 5.18. Let C be a link with number of components p > 2, then 
the Alexander polynomial A,c(l) = 0. 

To prove that the above polynomial is of Hosokawa type (denote it by Ac(t)) 
we use the following facts from knot theory [158]: if the Alexander polynomial 
is of Hosokawa type then, a) the multiplicity of the factor (t — 1) is even, b) 
Ac(t) is a Laurent polynomial of even degree (up to multiplication by factors 
t m with m — 0, ±1,±2, ...), c) the number of components p of the link should 
be greater or equal to two but otherwise is arbitrary Rcf. [157], page 121. 
In the present case Milnor [26] (and also A'Campo [156]) had shown that the 
number p of components in the link is given by Eq.(5.33)).; d) the Hosokawa 
type Alexander polynomial possess the property: Ac(t)=Ac(\)- Looking at the 
Alexander polynomial given by Eq.(5.35) it is straightforward to demonstrate 
that, indeed, the equivalence Ac(t)=Ac(j) holds. 

Remark 5.19. The above arguments apply, strictly speaking, only to 3d 
links embedded in S 3 . The Alexander polynomial obtained by A'Campo [156] 
describes, however, knots/links embedded not into S 2n+1 but into CP™. This 
circumstance creates no additional problems locally (since locally CP™ D C™ +1 D 
R 2n+2 D S 2 " +1 ). Nevertheless, the results related to the Alexander polynomial 
derived by A'Campo should be treated with some caution. We hope, that the 
Theorems 5.14. and 5.15. provide the necessary assurance that, indeed, it is 
possible to extend results for links in 3d to links of dimensions higher than 3. 

Remark 5.20. The Alexander polynomial, Eq.(5.32), formally looks differ- 
ent from that given by Eq.(3.20). This difference is superficial however since 
Eq.(5.32) is written for a special case of the Bricskorn-Pham singularity : when 
all a^s — N. In this case all roots of unity factors in Eq.(3.20) originate from 
the same Eq.(A.3) (with fixed n = N). 

Remark 5.21. The Corollary 3.7. should be amended now in view of 
Eq.(5.32). The mass spectrum is determined actually by two parameters: n 
and N. This option, is only plausible but, in fact, is less convenient than that 
discussed below. 



A(t) = (1 

= (1 
= (1 



tf[{t-p){t-- P )f 
tf(t 2 +t+ 1) 3 



(5.35) 
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5.3.2 Newton's polyhedra and the Hodge spectrum 



Based on the results of Appendix D, we conclude that the above two parameters 
control the Hodge spectrum of the B-P singularity. As results of Ref. [36] , page 
328, indicate, the total number of differential n+1 forms Lo b of the type 

6, b n+1 dsi ds n ds n+ i 
oj b = si 1 ■ ■ ■ s^+{ A • • • A A (5.36) 

is given by the Milnor number, Eq.(5.33), (e.g., please, read in addition the 
Remark 5. 26. below). Hence, this number provides as well the total number of 
possible Veneziano amplitudes. This number has a very interesting geometrical 
and topological interpretation which we would like to discuss now. We begin by 
introducing some definitions. 

Definition 5.22. The function x k = x^ 1 ■ ■ ■ x„™ is called monomial with 
the (vector) exponent k=(fci, k n ), ki <G Z + .The monomial is called Laurent 
ifheZ [158]. 

Let A C Z k be a finite set of integer vectors. Then, this set can be identified 
with the set of corresponding monomials. Let C" 4 be the space of Laurent 
polynomials made out of monomials from the set A, i.e. 

/(x) = f( Xl , ...,x n ) = «ka; k . (5.37) 

Evidently, the B-P variety, Eq.(3.8), is just a special case of the variety defined 

V(f) = {z e (C*) k | /(z) = 0} (5.38) 

where C* = C \ 0. Such variety is known in literature as the algebraic toric 
variety [158] .The name "toric" will be explained shortly. 

Definition 5.23. The Newton polytope (polyhedron) N(f) associated with 
function /(x) defined by Eq.(5.37) is the convex hull in R™ of the set {k : ^ 
0}. 

Remark 5.24. Use of the above definition of N(f ) in actual calculations 
may or may not be convenient. There are several algorithmically different (but 
combinatorially equivalent) ways to define a convex polytope in R™[159]. Ac- 
cording to Ref. [159], a convex polytope V is either 

a) a convex hull of finite set of points in R™ (which is essentially just the 
Definition 5.23) or 

b) is the result of an intersection of finitely many half spaces in R™. In the 

last 

case V is defined by the set if inequalities: 

P = {xeR": (U, x) > -Oi , i = 1, to} (5.39) 

The vectors U E (R™)* are from the dual space to be defined shortly below and 
a, G Rj, i = 1, to. Let H^ be an affinc hypcrplane in R™ then the intersection 
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"PlHHi is called face of the polytope. The boundary dV of V is the union of 
faces. 

Let us discuss now the rationale for introducing the dual space. To this 
purpose let us recall that the function / is considered to be quasi homogenous 
of degree d with exponents l\,...,l n if 

f(\ h x u \ ln x n ) = X d f{x u ...,x n ) (5.40) 

provided that A <E C*.The exponents k associated with variables called 
weights. 

In the light of the above definition, it is clear that each monomial in Eq.(5.37) 
is associated with the scalar product of the type: 

Y,(h)ikj = di (5.41) 

3 

with index i numbering the monomials. Hence, the vector 1 = (Zi, ...,l n ) is the 
dual of the vector k. The name "toric" comes from the observation that A 
can be chosen in such a way that |A| = 1. Hence, the usual topological torus 
T" (which is a product of n circles) defined by 

T" = {(e 2 ^,...,e 2m ^) G (C*)"}, 

with ip = (ip l7 ...,</?„) running through R", is just a subgroup of the algebraic 
torus. The algebraic torus (C*)™ is the semigroup with respect to component- 
wise multiplication [158], i.e. if t e (C*)" then, the action of t on x k is given 
by 

t 1 • x k = t 1 k x k 

which is just a special case of Eq.(5.40). Also, Eq.(3.26) is just a special case 
of this equation. The scalar product 1 • k = d defines a set C of hypcrplanes. 
Following Ref.[160] we call the set of all exponents of monomials that appear 
in /(x), Eq. (5.37), the support. It is clear that the members of the support lie 
either inside of the N(f) or at its faces. The hyperplane set C plays a major 
role in what follows. 

Remark 5.25. It can be shown [161] that the hyperplane set C determines 
convex polyhedron V completely and vice versa (i.e. they are dual to each 
other). Since the polyhedron N(f) defines(determines) toric variety, Eq.(5.38), 
this means that the hyperplane set C determines the toric variety as well so 
that all results of this work can be deduced from (encoded in) some 
properties of the hyperplane set C. The theory of these sets is known 
in the literature as the theory of arrangements. Condensed, physically relevant 
and pedagogically clear exposition of this theory can be found in the book by 
Orlik and Terrao [162]. This book among many other things discusses work of 
Schchtman and Varchenko [137] connecting theory of hyperplane arrangements 
with the theory of K-Z equations 23 . 

23 The systematic exposition of K-Z equations involving Orlik-Solomon algebra related to 
hyperplane arrangements can be found in Ref.[78]. Alternative treatment can be found in 
recent MIT lecture notes by Varchenko [80] . 
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With such background we are ready now to go back to the differential n+1 
form, Eq.(5.36), in order to rewrite it in the language of toric varieties. It should 
be clear that u>b in Eq.(5.36) can be presented as follows 



UJ b = x b w 



(5.42a) 



where b=(6i, b n +i) and ujq is given by 




dx n+ i 



X n +1 



(5.42b) 



Taking into account Eq.s(3.26) and (5.41) and using Eq. (5.42a) we obtain the 
hyperplane equation which, in view of Eq. (5.20b), can be written as 



and should be compared with Eq.(5.39). 

In writing this equation we have taken into account that the eigenvalue A of 
the monodromy operator (Section 3.2) is related to a as 



in accord with notations of the Appendix D. 

Remark 5.26. From Eq. (5.20a) it should be clear that the vector b is just 
one of many possible allowed vectors. The totality of these vectors is given 
by the Milnor number fj,, Eq.(5.33). The set of these vectors is in one-to-one 
correspondence with the set of eigenvalues {A} of the monodromy operator, 



For the B-P polynomial, Eq.(3.7), it is easy to give useful geometric descrip- 
tion of these vectors. Indeed, consider a hypercube (0,N) n+1 = (0, N) x ■ ■ ■ x 
(0, N). Then a vector b is represented by one of the points strictly inside of the 
hypercube. If the volume of the unit cube cell is chosen to be unity, then the 
number of points inside the cube coincides with the volume (N — l) n+1 which 
is just the Milnor number /i,Eq.(5.33). 

In view of the results of Appendix D, the totality of hyperplanes described by 
Eq.(5.43) is in one-to-one correspondence with the totality of the Hodge numbers 
associated with the differential forms, Eq.(5.36).We would like to describe this 
connection in some detail now following the work of Arnol'd [163]. 

From Eq.(D.6) it follows, that -a in Eq.(5.44) should be a positive rational 
number. Mathematically, both integer and non integer numbers are allowed. In 
the case of integers the corresponding monodromy eigenvalue A — 1 , while in the 
non integer case A ^ 1. In view of Eq.(3.17a), only non integers should be of 
physical interest. Physical relevance of the integer case requires further study. 
So, let -a be some integer + a fraction (in Appendix D it is denoted as a). The 
following number-theoretic problem can be posed now: For the set of integers b 




(5.43) 




(5.44) 



Eq.(3.17). 
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lying inside the hypercube (0, iV)™ +1 and for positive rational number a lying 
in the range 

n — p < a < n — p+1 (5.45) 

(with p being defined either by Eq.(5.9) or in the Appendix D) find how many 
vectors b produce -a in Eq.(5.44) provided that a lies in the above domain. 

The Hodge number h pq can be found now if in addition to the number- 
theoretic problem just described we would find the number of solutions of 
Eq.(5.43) for case when 

—a = a = n — p+1. (5.46) 

Denote h 1 ^ 1 the number of solutions for the first problem and h pq the number 
of solutions for the second problem then, the Hodge number h pq is given by 

h pq = hf + h\\ (5.47) 

e.g. see Eq.(D.lO) of Appendix D. 

Geometrically, h p ^ is just the number of integer points in the hypercube lying 
between the hyperplanes defined by the left and right sides of the inequality, 
Eq.(5.45), while h\ q is the number of integer points lying in the hypcrplanc 
defined by Eq.(5.46). According to the Ref.s[159] and [163], N(f) for the B-P 
polynomial is combinatorially equivalent to the hypercube (0, N) n+1 = A whose 
vertices belong to the integral lattice Z™ +1 C R n+1 . For any positive integer k, 
let l(kA) be the number of integral points in fcAnZ™ +1 . By definition, 1(0 A) = 1. 
In addition, let l*(kA) be the number of integral points in the interior of fcA. 
Based on these definitions, two generating functions (Poincare' polynomials) 
can be constructed now: 

P A {t) = J2 l (k&)t k (5.48a) 

k>0 

and 

Q A {t) = J2l*(kA)t k . (5.48b) 

fe>0 

Batyrev [164] proved that, if additionally one defines two related functions : 

*a(*) = (1 -t) n+2 P A (t) = ]T^(A)f (5.49a) 

and 

*a(*) = (1 - t) n+2 Q A (t) = J2 Vi{*)t\ (5.49b) 

i>0 

then, 

t n+2 *A(* -1 ) = *a(*). (5.50) 
According to Arnol'd, Rcf. [163 ], for the B-P variety one obtains 

/ i \n+l 

Qa( * )= (it^) • (5 - 51) 
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Consider now the limit: t — > 1 +e, e —* + ,of the above generating function. An 
easy calculation produces: 

Qa(1) =M= (A^-l) n+1 , (5.52) 

where, again, it is Milnor number, Eq.(5.33). According to Eq.(5.33) this result 
also allows calculation of the Euler characteristic x{&n) of the complement 
of the Fermat hypersurfacc in CP™and, hence, the Euler characteristic of the 
Format hypersurface as well [156]. 

Remark 5.27. The result, Eq.(5.52), is known in the literature as Koush- 
nirenko theorem [165]. 

Remark 5.28. Based on the observation that [166] for any Kahler manifold 

X 

n+l 

X (X) = {-l) v+q h p \X). (5.53) 

p,q=0 

Danilov and Khovanskii [167] had proposed to look for properties of the gener- 
ating function 

e(X; x,x) = Yj e pq (X)x p x q , (5.54) 
p,q 

where e pq (X) = {-l) p+q h pq (X). Because of the connection between N(f) and 
X exhibited in Eqs.(5.33) and (5.52), these authors had demonstrated how to 
recover the Hodge numbers h pq from generating functions, Eq.(5.48), of the New- 
ton polyhedron thereby connecting the generating function e(X; x, x) with that 
for the Newton polyhedron. Useful condensed summary of results of Rcf. [167] 
can be found in the review paper by Cox [168] . 

Remark 5.29. The formula of Arnol'd for Qa(£), Eq. (5.51), is a special 
case of more general formula obtained by Steenbrink [169], Brieskorn (unpub- 
lished) and Saito [170]. Their results are alternative to that obtained by Danilov 
and Khovanskii and, in our opinion, more efficient computationally. In the case 
if /(x), Eq.(5.37), is linear combination of monomials x k such that Eq.(5.41) 
acquires the form lk = 1 for all monomials entering /(x), it is possible to 
demonstrate that 

11+1 t li — t 

Ga(*) = 11(1— T> ] (5 ' 55) 

i=i 

which for the case of B-P polynomial is reduced back to Eq. (5. 51). However, the 
interpretation of the polynomial Q&{t) by Saito [170] is different. According 
to Saito 

Qa(*) = J> Q< 

i 

wither = — InAi (e.g. see the Appendix D and Refs. [135,170]). Hence, Q^{t) 
is also a generating function for the Hodge spectrum. Since earlier (and in the 
Appendix D) we had mentioned that the spectrum is degenerate, the above 
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generating function acquires the following form: 



{a} 



.a 



(5.56) 



with degeneracies n a = J2 n a,i- The numbers n Qi ; participate in yet another 



with the spectral pair (a, I) being defined in the Appendix D. Since both a and 
/ in this pair are related to p and q, one can identify n a j with h 1 ^ 1 [135]. Hence, 
the generating function Qa(£) allows, in principle, determination of the Hodge 
numbers h p ^ . Examples are provided in Ref.[135]. The same results (with much 
less details) can be found in a short paper by Varchenko and Khovanskii [171]. 



6 Discussion 

To our knowledge, so far the results from geometry of toric varieties had been 
used in the so called mirror symmetry claculations [57]. The results presented 
above indicate that actually all results of conformal field theories and high en- 
ergy physics could be described from the standpoint of geometry of toric vari- 
eties. This discipline itself is not new however. Actually, it can be considered 
as branch of the theory of the linear algebraic groups to which some results of 
the theory of Lie groups and Lie algebras had been adopted. For instance, the 
result, Eq.(5.56), has its analogue in the theory of characters in the theory of 
Lie groups/ algebras and, hence, one can derive the Weil formula for characters, 
etc. Under such interpretation, the Milnor number /i, Eq.(5.33), plays a role 
of the dimension of representation of the Weil/Coxeter reflection group. With 
little extra effort one can reobtain Kac-Weil character formula-cental for de- 
velopments of conformal field theories. In addition, obtained results allow the 
symplectic interpretation enabling us to restore the underlying physical model 
producing , for instance, the Veneziano amplitudes. The results of Atiyah doc- 
umented in Ref. [40] indicate that this model is also going to be a string. The 
combinatorial aspects of the Lie group theory representations will allow us to 
restore the Verlinde fusion rules as a special case of Schubert calculus consid- 
ered in our earlier work on Kontsevich-Wittcn model [64]. All results of this 
earlier work and the work just presented can be reobtaincd with help of the 
symplectic formalism based on the Duistermaat-Heckman formula [40]. If one 
is interested in development of p-adic aspects of string theory, the concept of 
motives and its latest development culminated in the works by the latest Field 
medalist, Vladimir Voevodsky, can be found in the form accessible to physicists 
in Refs. [186,187]. Details of the results outlined above are to be presented in 
subsequent publications. 



generating function 




(5.57) 
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Appendix A. Some results from cyclotomic field theory 

Suppose we have the number field F p of characteristic p, that is the coset 
Z/pZ . Let Ci € F p . It can take p values. Suppose now that we would like 
to extend the field F p . This can be accomplished very easy using concepts of 
standard linear algebra. To construct the vector space of dimension n we need 
to have a set of n linearly independent (basis) vectors ej so that any vector A 
can be decomposed as usual: 

n 
i=l 

Let now the role of <2j is to be played by Cj and the role of -by some basis 
elements on to be determined momentarily. Then, by construction, a number 

n 
i=l 

belongs to the extension of the field F p . Since each a can have n values, the 
field F p [a] contains p n elements and is known in literature as the Galois field 
GF(p").The task remains to find the basis set explicitly. The familiar example 
of complex numbers helps a lot in this respect. Indeed, the complex number 
i = \/ —1 is just solution of the equation x 2 + 1 = 0. Clearly, all complex 
numbers z are given by z = al + ib with a and b belonging to R, Q or Z. The 
basis elements are 1 and i. To go beyond this elementary case it is sufficient to 
consider all solutions of the equation 

x n -x = 0. (A.3) 

It can be shown [99] that the above equation can be also rewritten as 

n 

x n - x = Y[(x - a,) (A.4) 

»=i 

with on being a root of Eq.(A.3). In addition, it can be shown that all roots 
are distinct and, hence, there are n of them. Naturally, they can form the basis 
which we are looking for. Moreover, Eq.(A.3) can be equivalently rewritten as 

x™- 1 = 1 (A.5) 

and, hence, the cyclotomic (that is circular) nature of this equation becomes 
clear. Indeed, consider Q k — exp{ia} with a = k-^j. If for k we assign values 
0, 1, 2, n — 2, we obtain n — 1 equidistant points on the circle. It is clear that 
this construction can be performed for any n and, hence, if we consider instead 
of Eq.(A.5) the equivalent equation 

x n = 1 (A.5a) 
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we shall obtain exactly n values for £ fc . These are independent of each other 
and can be used instead of a-s in Eq.(A.2). Moreover, there actually no need of 
substituting all Q's into Eq.(A.2). It is sufficient to use just one of the so called 
primitive roots of unity in Eq.(A.2). By definition all roots of Eq.(A.5a) can 
be obtained as powers of those which are considered to be primitive. To obtain 
primitive roots the following steps should be taken. For a given n one should 
look for all k's such that k \ n and k < n . The number of such possibilities is 
given by the Eulcr's function (p(n). Recall that the coset Z/nZ is a field (and is 
denoted F„) only if n is some prime number p. An element Ci is a unit in F„ 
if (ci , n) = 1 (that is a \ n). Hence, there is one-to-one correspondence between 
the units in F„ and the primitive roots in cyclotomic fields. So if n = p we have 
<p(p) = P— 1 primitive roots. If we take one of them, say C, then the desired basis 
is given by :1, C, C 2 > C P ~ ■ Let us renumber these primitive roots as follows : 
Ci) C2) Cip(p)- Moreover, let us remove for a moment the restriction on n to be 
a prime, that is we assume that n=n"=i r».Then, surely, <^(n) = <p(ri) ■ ■ ■ f(r m ) 
and for each <p(ri) the cyclotomic polynomial <J> ri (x) can be defined now as 

(X-Cl)---(X-C vi r i )) = ^r i {x). (A.6) 

It is monic irreducible polynomial of degree <p(ri) by construction. Using multi- 
plicative property of Euler's function it is clear that the following identity should 
hold 

x n -l = Y[$ d (x) (A.7) 

d\n 

with d running over all divisors of n. For n — p the result can be considerably 
simplified by noticing that x = 1 is the root of Eq. (A. 7). This means, in view of 
Eq.(A.6), that Eq.(A.7) can be simplified to 

x p - 1 = (x - l)^- 1 + x p - 2 + ■ ■ ■ + x + 1). (A.8) 

By construction, the second multiplier in the r.h.s of Eq.(A.8) is monic irre- 
ducible (in R, Q or Z) polynomial. It should be clear that the roots of the 
equation 

x p - 1 +x p - 2 + ---+x + l = (A.9) 

are primitive roots of the cyclotomic Eq.(A.5a). Thus, the connection between 
the cyclotomic fields and polynomials is clear. This connection can be extended 
further with help of Eq.(A.2). Indeed, let us write 

n-l 
i=0 

and also 

i=0 
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etc. and, finally, 

71-1 

i=0 

for some cti.j in F p . Now we want to look at equation 

n 

^ Xl (i l =0 (A.10) 

i=0 

with Xi GF p . Upon substitution of equations for powers of (i' s in Eq.(A.lO) 
the representation of rr-s through a-j s can be found and, hence, to construct 
equation of order n with j3 being a root. The polynomial associated with such 
an equation is irreducible and always can be made monic [49], page 430. Since 
for a given prime p there are p n (3's according to Eq.(A.2), hence, there are p n 
monic irreducible polynomials of degree n. If we exclude the trivial case (3 = 0, 
then there are p n — 1 nontrivial irreducible polynomials. 

Remark A.l. Let S m be symmetric group acting on complex space C m 
by exchanging coordinates. Then 

Y m := {(01, z m ) eC m \ Zi ^ Zj for i ^ j}/S m 

is the set of all unordered m-tuples of distinct complex numbers. It can be 
shown (Bricsk, malle &matzat) that the set Y m is in one-to-one correspondence 
with the set of monic irreducible polynomials with roots z\,...,z m . Define a 
braid on m strings (and initial point z°, z°J as homotopy class of a closed 
paths in Y m with initial and final point at z®, z^ . Such constructed braid is 
in one-to-one correspondence with the set of monic irreducible polynomials 24 . 

Appendix B. Some results from Nielsen-Thurston theory of surface 
homeomorphisms 

Let S be closed orientable Riemann surface of genus g.The first homotopy 
group, the fundamental group 7Ti(iS) of surface S is made of 2g generators 
{xi, yi}, i = 1 — g and a single relation so that its presentation is known to be 

tti(5) =< xi,yi, :.,x g ,y g \ [a^yi] • • • [x g ,y g ] > . (B.l) 

Nielsen has noted that there is one to one correspondence between automor- 
phisms of 7Ti(5) and surface self-homeomorphisms. This is summarized in the 
following proposition 

Proposition B.l. (Nielsen [172]) If g > 1, then every element of Out(iTi(S)) 
is represented by a unique isotopy class of self-homeomorphisms of S. 

24 Braids arc connected directly with the hyperplane arrangements mentioned in Section 
5.3.2. in connection with the K-Z equations. 
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An important subgroup of Out(iri(S)) is the mapping class group M g . 
Geometrically, this group is finitely generated by the Dehn twists in simple 
closed curves (lamination set) on S whose physical significance was discussed 
extensively in our previous work, Refs[4-6]. A simple closed curve C on the 
orientable surface S has a neighborhood £ homeomorphic to an annulus which 
is convenient to parametrize by {[r,#]|l< r < 2}. The Dehn twist in C can be 
imagined as an automorphism Tc : S — > S. It is given by the identity off £ and 
by ► [r,6 + 2irr] on £ . We would like to illustrate these concepts on the 

simplest example of the punctured torus T 2 using results of our previous work, 
Ref. [4]. In this case Out(m(T 2 )) = GL 2 (Z) and Mi,i = PSL{2, Z). Since any 
transformation from PSL(2, Z) is obtainable by projectivisation of SL(2, Z) we 
discuss everything in terms of SL(2, Z) with projectivisation at the end. Any 
transformation which belongs to SL(2, Z) is expressible in terms of 2 x 2 matrix 
A given by 



with integer coefficients subject to condition:detA=afe — cd=l. The character- 
istic polynomial for this matrix is given by 



This implies that the eigenvalues of A are either : 

a) both complex ( when trA=0,l,-l), 

b) both equal to ±1 (when trA=±2), 

c) distinct and real (when \trA\ >2). 

If $a is toral automorphism then, transformation a) is called periodic since 
(3t0" = 1 f° r some n (actually, n = 12 by the Hamilton-Cayley theorem) 
, transformation b) is called reducible since it leaves a simple closed curve C 
invariant, transformation c) is called {pseudo) Anosov if the (line) vector field 
on surface (does) does not contain singularities. Physical significance of this 
difference is explained and illustrated in our previous work, Refs.[5,6]. 

The largest of two eigenvalues is associated with the topological entropy of 
the (line)vector flow and is related to the amount of stretching of surface S and, 
hence, with the dilatation parameter of the Teichmullcr theory. The Nielsen- 
Thurston theory generalizes the above classification of surface automorphisms 
to all surfaces of genus g > 1. Already Nielsen had realized [172] that for g > 1 
it is more convenient to study homcomorphisms of surface S by considering 
their image on the universal cover of S which we choose as Poincare' disc model 
of H 2 , i.e. int DUS^ =H 2 . According to Nielsen [172] 

Proposition B.2. Any lift h of the surface self- homeomorphism h: S — > S 
to the universal cover of S extends to a unique self-homeomorphism of the unit 
disc D, i.e. to int D U S^. 

Surface sclf-homeomorphisms h are associated with the Dehn twists con- 
nected with a set of simple closed nonintersecting curves homotopic to geodesies 




(B.2) 



t 2 - trA t + det A = 0. 



(B.3) 
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(such set is called geodesic lamination C). Their lifts h(C) are associated with 
some maps of the circle extendable (quasi conformally) to the interior of the 
disc D e.g. see Ref.[73]. An image of the closed geodesic on S, when lifted to 
H 2 , is just a segment of a circle whose both ends lie on S^. Since geodesies are 
nonintersecting, circle segments on are also nonintersecting. In this work 
we are interested only in the periodic maps of the circle since non periodic maps 
can be always systematically approximated by periodic ones as it is explained 
in terms of continuous fractions and associated with them Dehn twists in our 
earlier work, Ref.[4] or, alternatively, in Milnor's lecture notes [173] and Ref.[55]. 
In connection with such circle maps the following remark is of importance. 

Remark B.3. (A variant of Sarkovskii theorem, Ref.[174], page 88). Let 
/: S 1 — ► S 1 be a continuous map of the circle with a periodic orbit of period 3. 
If the lift /: R — > R has also a periodic orbit of period 3 then, / has periodic 
orbits of every period. The condition on the lift map f cannot be dropped. The 
continuous maps of the circle can be replaced by the piecevise linear maps while 
3 still remains as minimal period. 

Remark B.4. As noted by Kontsevich [175], the moduli space problem 
makes sense only for Riemann surfaces obeying the following set of inequalities 

g > 0, n > 0, 2 - 2g - n < (B.4) 

with n being the number of distinct marked points (effectively, distinct bound- 
ary components). Boundary components can be eliminated by the Schottky 
double construction. This construction can be performed as follows. If M is 
a complex manifold with Ci, C n boundary components, one can consider an 
exact duplicate of it, say M, with the same number of boundary components, 
say, C\,...,C n . Evidently, for each point x € M there is a "symmetric" point 
x e M . The Schottky double 2M is formed as a disjoint union MUM and iden- 
tifying each point x € Ci with point x e Ci for 1< i < n. In the simplest case 
we have initially cither punctured torus, i.e.g = 1, n = 1, or the trice punctured 
sphere, i.e. g = 0, n = 3. In both cases the Schottky double is a double torus. 
A double torus has 3 geodesies which belong to the geodesic lamination C. The 
image of these geodesies lifted to H 2 produces 3 circular arcs whose ends lie on 

. This is the minimal number of arcs required for the moduli space problem 
to make sense. According to Remark B.3. this is also the minimal period for 
the periodic homeomorphisms of the circle in view of the Sarkovskii theorem. 

In Ref.[176] it is argued that the total number of geodesies on the Schottky 
double is 6g — 6 + 3n. This is the dimension of space of holomorphic quadratic 
differentials (real on each of the boundary components). Hence, in accordance 
with Teichmiiller theory [8] , it is the dimension of the Teichmullcr and, accord- 
ingly, the moduli space of such Schottky doubled surface. 

Appendix C. Complex multiplication and Hecke operators 

In Section 4b we have introduced both lattice V and sublattice L . This 
caused us to introduce the relationship: ad — be = n,n > l.We would like now 
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to study physical implications of this relationship. To this purpose, following 
Rcf.[89], we introduce 3 types of matrices: 



T> n = " b d ) £ M-AZ): ad - be-, n 



S n = b d ^j eM 2 (Z); ad = n,a,d> 0,0 < 6 < dj 

° J jeSL 2 (Z); a,/5, 7 ,<5€M-/5 7 =l 

Naturally, S n < — > SL2CZ) \ T> n . Introduce now the homothety operator R\ as 
follows 

R X L = \L (C.l) 

where A e C and the action of n th Hecke operator T(n) on the lattice L is 
defined by 

T(n)L= i L ')- (C-2) 

lcl' 

[L:I/]=n 

In order to use this equation, it should be rewritten in more explicit form as 

T(n)L= ]T a(L). (C.3) 

aes n 

It can be demonstrated [89] that 

a) R\R^ = R Xfl e C; 

b) R\T(n) = T(n)R x VA e C,n > 1; 

c) T(mn) = T(m)T(n) Vm, n > 1 with gcd(m, n) = 1; 

d) T{p e )T(p) = T(p e+1 ) +pT(p e - 1 )R p for p prime, e > 1 
and p \ n , while for p \ n we get 

d')T(p*) = [T(p)] e . 

Since any number n can be decomposed into primes the above rules are 
sufficient for any n, prime or not. From the definition of the Hecke operator 
it is clear that it is acting on the lattices, more exactly, it acts on equivalence 
classes between different lattices. Let now f(L) be some function defined on 
such lattices. We shall assume that this function is homogenous of degree -k. 
That is 

/(ZWI + ZC0 2 ) = LU^fiLUl/u,). (C.4) 

The requirement of homogeneity makes this function invariant with respect to 
changes of the basis. Indeed, if we have change of the basis given by 



? - 



a b 

c d J \ u)i 
then, 

/(Za>; + Zw' 2 ) = w'fV^aM) = (owi + doj 2 )- k (d + c^/uj^f /(w 2 /wi) 
= ^fi^/uj,), (C.5) 
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provided that / is modular function of weight k, that is 

/(( ! \ )r) = (cr + d) k f(T) for ( " \ )eSL(2,Z). (C.6) 



c d J ' ^cd 

As consequence of Eqs.(C.4) and (C.5), function f(L) possess the homogeneity 
property 

f(aL)=a- k f(L). (C.7) 

Consider now action of the Hccke operator on function f(L). In the case if such 
a function is a modular form of weight k we obtain 

(T fe (n)/)(r)=n fe - 1 ]T f(L')=n k ~ 1 £ d~ k f(^), (C.8) 

j r j' ad—n,a>l 
[L:£']=n "< b < d 

where, following Ref.[89], we defined (T k (n)f) (r) = (n* -1 !^)/) (L). The 

result, Eq.(C.8), was obtained with help of Eqs.(C.3)and.(C6). In actual cal- 
culations one should take into account that the weight k should be only even 
number. The modular forms with odd weight are zero as can be easily shown. 
The factor n k ~ 1 looks somewhat artificial at this point. Hence, we would like 
to explain its origin now. 

The modular function /(r) admits the Fourier decomposition /(r) := c(m)q m , 
where q — exp(27rir). The eigenvalue problem 

(T k (n)f)(r) = X(n)f(r) (C.9) 

in view of Eq.(C.8) can be written as follows 

A(n)/(r) = n k - x d ~ k J2 c ( m ) cx p( 27 "™K + b)/d) (CIO) 

ad— n,a>l mEZ 
0<6<d 

= n*" 1 Yl c(m)d- k e 2 ™ maT / d e 2mmb ' d . 

m£Z ad=n,a>l 0<b<d 



Since 

d if d | m 
if d j to 



^2nimb/d 



E 

0<b<d 

we have to replace to by md = ran I a thus producing n/d — a resulting in 
cancellation of the prefactor n fe_1 . After this, Eq.(C.lO) acquires the following 
form 

A(n)/(r) = ^ ^ a fe-i c (!!^) e 2™»ar^ (cn) 

m G Z ad— n,a>l 

The last result can be further rearranged thus producing 

A(n)c(/)= £ « fe ~ lc (^)- (C- 12 ) 

a|gcd(/,n) 
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Suppose gcd(7,n) = l,then, surely, a = 1 which leads us to X(n)c(l) = c(nl). 
Let further, I = 1. This produces :A(n)c(l) = c(n). From here it follows that 
c(l) 7^ and it can actually be taken as 1, Ref. [89], page 78 and Ref.[100], 
page 101. If this is so, then we are left with a very nice result: X(n) — c(n), to 
be used in the main text. Moreover, we will also need the following result: 

(Rj) (L) = f(XL) = \- k f(L). (C.13) 

Using this result in equation 

(T(p e )T( P ).f) (L) = (T(p^)f) (L) + (pT(p^)Rj) (L) (C.14) 

we obtain 

(T(p e )T(p)f) (L) = (T(p e+1 )/) (L) + (p 1 -*^- 1 )/) (L). (C.15) 

Using the definition {T k {n)f) (r) = (n k - 1 T(n)fj (L) in Eq.(C15) and multi- 
plying both sides by factor p( e+1 )( fe_1 ) produces our final result 

T( P e )T(p)f = T(p e+1 )f + p'-'Tipe- 1 )/ (C.16) 

to be discussed from yet another (physical) point of view immediately below. 

To this purpose we would like to remind our readers some excerpts from the 
graph theory. A graph G consists of finite set of vertices, V and edges E along 
with the following two maps: 

E xV ande^ (o(e),t(e)) (C.17) 

so that Ve <G E we may have e — > e which means that the edge can change the 
orientation into opposite so that o(e) = t(e). The element o(e) is the origin of e 
and t{e) is the terminus of e.Two vertices Vq and Vi are adjacent if there is an 
edge with o(e) = i>o and t(e) = Ui.The order of G, denoted as \G\ ,is determined 
by the number of vertices in G. A walk W in G is a sequence of oriented edges 
ei, e r such that i(ej) = o(ej+i) for 1 < i < r — 1. If walk is from t>o = o(ei) 
to v r = t(e r ) it means that the walk is from vq to v r of length r. We shall be 
interested only in walks without backtracking . For such walks, let \G\ = n 

be the order of the graph then, the matrix A r =(af[J^ is a symmetric n x n 
matrix with nonnegative integer entries counting the number of walks of length 
r from Vi to vj. By definition, A\ is the adjacency matrix [177]. The degree 
degVi of the vertex Vi is the number of edges emanating at Uj. Let D be the 
diagonal matrix with entries deg vi, deg n n along the diagonal then, matrices 
A r can be determined recursively from the following set of equations 

A X A X = A 2 + D, (C.18) 
A r Ai = A r+ i + A r -i(D - I) for r > 2. 
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The proof of this result can be found in Ref.[178]. In order to use these recursions 
we would like to make a few simplifications. First, we would like to consider 
only fc-regular graphs (i.e. such for which degree of all vertices is the same and 
equal k). Next, let k = p + 1 for some prime p. And, finally, let A r = B r — £> r _2 
for all r > 1. Then, the system of recursions given by Eq.(C.18) is reduced to 



provided that B-\ = 0, Bq = I, B\ = Ai.This result should be compared with 
Eq.(C.16) (for k = 2). Evidently, the total agreement is reached. It should 
be clear that the edges of G represent equivalence classes between different 
lattices represented by vertices of G. The universal covering of such graph is 
p + 1 -valent tree known in the literature as Bruhat-Tits tree [179]. Should 
our lattices come from the lattices whose base is made of more than 2 elements 
then, instead of trees, we would have buildings [179, 180]. Such necessity indeed 
occurs for higher dimensional complex tori (e.g. read Section5) and had been 
discussed in connection with geometry of toric varieties (Section 5) by Mumford 
[181] already in 1970. 

Because of the agreement between Eqs.(C.19) and (C.16), it is sufficient, 
in principle, to do all relevant physics "locally", that is for the spe- 
cific prime p, in order to restore the "global" picture for all primes. 
Although this idea is not new, as discussed in the Introduction, its implementa- 
tion in the present context may still require considerable efforts. E.g. to rewrite 
this paper using the terminology of motives [38,145] would require entirely new 
(and much more advanced) level of presentation. This task is left for the future. 

For now we only would like to notice that since A\ is the adjacency matrix, 
this means that the Hecke operators can be identified with the adjacency matrix 
of some graph G.This observation makes them already physical especially be- 
cause the discrete (combinatorial) Laplacian A for /c-regular graph G is known 
to be 



Harmonic analysis for such type of Laplacians is discussed in detail in Refs. [16,17,182- 
184]. Higher dimensional analogs of Hecke operators, Hecke algebras, etc. are 
discussed in Ref.[185]. 

Appendix D. Veneziano amplitudes from the oscillatory integrals. 

In this appendix we would like to provide the condensed summary of known 
results for oscillatory integrals [36,144]. Fortunately, these results are sufficient 
for reobtaing the Veneziano amplitudes discussed in the main text. 

Following Ref.[36], Ch-r 11, we are interested in study of integrals of the 
type 



B r B\ = B r+ i +pB r -i for r > 1 



(C.19) 



A = M-A 1 . 



(C.20) 




(D.l) 
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where the function /: (C n ,0) — > (C,0) is holomorphic in a neighborhood of 
its only one critical point located at the origin, [r] is n— dimensional homology 
chain associated with the cohomological differential n— form w. The parameter 
r is expected to be positive. At the boundary of the chain [r] the real part of the 
function f(z) is negative. The chain [r] with such property is called admissible. 
We shall need the following identity (Lemma 11.2. of Ref.[36]) valid for large 
t's : 

to f \ 

e Tf Wu> w J e- Tt J co/df dt (D.2) 

[r] o \ a _ t[ r] / 

with to being a some positive number which eventually can be taken to the 
infinity. This result can be easily understood in physical terms. Indeed, let us 
consider an auxiliary integral 

I{t)= J 5{t-f{z))w (D.3) 
[r] 

with 6 being the usual Dirac's delta function. Such type of integrals formally 
occur in many calculations, e.g. in many- body quantum mechanical problems 
involving the density of states calculations, especially for the phonon spectra, 
etc. As it is shown in Ref.[148], Ch-r 4, 



cu/df (D.4) 

so that for t — > oo the r.h.s. of the Eq.(D.2) is just the Laplace-like transform 
of the density of states. 

Let now f(z) = z N+1 and u) = z l ~ x dz with 1 < I < N. Then, J(i) is given by 
(equation (7) on page 305 of Ref. [36] ) 

I(t) = const m^^ 1 (D.5) 

with const— y 3> . The phase factor Cj plays the same role here as in 
Sections 3.2. and 5. 2. Using this result in Eq.(D.2) and letting to to become 
infinity produces 

e Tf{z) uj w constT{-a)T a (D.6) 



[r] 

with a = j^pj- The result just obtained can be broadly generalized with help 
of the Fubini's theorem which can be stated as follows: 



/ 



[ri]x[r 2 ] [r x \ [r 2 
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If we combine this formula with the asymptotic result, Eq.(D.2), then, the 
Fubini's theorem becomes a statement from the theory of Laplace transforms: 
the Laplace transform of the convolution is equal to the product of the Laplace 
transforms. Looking at Eq.(D.5) it is easy to write a convolution. With accuracy 
up to a constant we obtain 

* 

I (a + P) = J ds\t- s\ & \sf = t &+f)+1 B{a + 1 ,p + 1) (D.8) 
o 

where a = — 1, /3 = — 1, 1 < I, m < N. Evidently, use of this result 
in the l.h.s. of Fubini's relation superimposed with Eq.s(D.2-D.5) produces 
back Eq.(5.24a) of the main text. Since the results just described and many 
others can be found in the comprehensive paper by Varchenko, Ref. [144], there 
is no need to repeat here his arguments related to Hodge and mixed Hodge 
structures associated with such type of oscillatory integrals. It also should 
be clear from reading of the same reference that the result, Eq.(5.32), for the 
Alexander polynomials is connected with Fubini's theorem just discussed. 

Following Refs. [160,163] we shall call the exponent a in Eq.(D.6) an order 
of the differential form u>. It is also called the spectrum of the critical point of 
the germ f(z) = z N+1 . Let now we have two germs / and g with respective 
spectrum {cti} and {Pj} so that i — 1, and j — l,...,rj. Then the spec- 
trum of the critical point of the germ / + g is obtained {cti + (3j + 1}. For the 
composition of germs leading to the B-P variety, e.g. see Eq.s (3.7)-(3.8), the 
spectrum is made essentially of numbers which belong to the set X(S 1 ) (actu- 
ally we have to divide these numbers by N) defined by Eq. (5.22a). In view of 
the Eq.s(3.17a) and (3.26) each number a is also associated with the eigenvalue 
A of the monodromy operator of the critical point, i.e. A = exp(±27ria). 

The dimensions of spaces of (p, q) differential forms lo defined by Eq.(5.9) 
are given by the Hodge numbers h p,q . These numbers had been related to the 
spectrum by Steenbink [169]. We follow, nevertheless, more recent and more 
clear exposition of this subject presented in Ref. [135] ( Ref. [144] is also very 
helpful). Let p + q = m,m < n, and choose a pair of numbers (a, I) where a is a 
spectral number a = — ^ In A normalized by the level p,i.e. n—p— 1 < a < n—p, 
and I is the weight number given by 

p + q — 1 ii A = 1 

Let h p ^ q be the dimension (in physics terminology, the degeneracy) of the 
eigenspace related to the eigenvalue A. Then, the Hodge numbers h p ' q are de- 
termined simply by 

h pq ^^2 h \ q - ( D - 10 ) 

A 

Taking into account the symmetry of the Hodge numbers (and also the numbers 



h™ = hl P ; hl q = h n-l-1,n-l- P {QlX ^ 1 ( D . lla ) 
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and 

h p,q = h n-q.n-p f Qr \ = 1 (D.llb) 

and using this symmetry the following "sum rule" is obtained: 

P,Q P,Q 

where the Milnor number \x is defined in the main text (for the special case of 
B-P variety only!) by Eq.(5.54). In addition, there is yet another sum rule 

5>=M^-1)- (D.12b) 

i 

Consider now the simplest example based on analysis of the germ f(z) = z N+1 . 
According to Eq.(D.6) we have a = j=p^ with 1 < I < iV.Clearly, \a\ < 1 
so that A^l,a= ivTT- ^ n y i cw °f the results presented above, we conclude 
that p = and , due to symmetry, also q = . This produces — 1 and 
h a >° = N = /i. Let us consider less trivial example of germ of the Fermat-like 
curve. It is relevant to the discussion in Section 3.2. In particular, we had 
obtained the set of differential forms, Eq.(3.24), distinguished by the numbers r 
and s (or a = and b= jj). Replacing now N by N+ land using Eq.(D.8) and 
the discussion which follows this equation, we obtain the spectrum as a + b — 1 . 
Clearly, the factor of -1 can be discarded since it does not carry any additional 
useful physical information. Hence, we are left with the spectral set with 
1 < m, n < iV.Using the arguments presented above the Hodge spectrum can 
be determined rather easily on a case by case basis. Hence, the particle 
spectrum is in one to one correspondence with the Hodge numbers 
of the spectral set. The alternative way of computation of Hodge numbers 
is presented in Section 5.3.2. 

We would like to conclude this appendix by writing down the set of the 
Picard -Fuchs equations (e.g. see page 336 of Ref.[36]) for the B-P type of 
singularity (Sections 3 and 5). They are given below: 

jjk jk 

-M=(r k -D T , (D.13) 

where the subscript k = (ko, k n ) numbers different integrals of the type given 
by Eq.(D.3) and is given by 

k + ... + k n 

Tk = 



N + l 

with 1 < k 0} k n < N. 
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